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, ^ I Abstract 

' The stability question of the Lane-Emden stationary gaseous star configurations is 

an interesting problem arising in astrophysics. We establish both linear and nonlinear 
dynamical instability results for the Lane-Emden solutions in the framework of the 

■ Navier-Stokes-Poisson system with adiabatic exponent 6/5 < 7 < 4/3. 

^ ■ 

■ 1 Introduction and formulation 

One of the simplest fundamental hydrodynamical models to describe the motion of self- 

^ ■ gravitating viscous gaseous stars is the compressible Navier-Stokes-Poisson system, which 

1> • can be written in Eulerian coordinates as follows: 
00 



in 



dtQ + div(^u) = 

9t(^u) -Fdiv(^u(8)u) +divS' = -pV$ (1.1) 



where (x, t) G x R"*", f)(x, t) > is the density, u(x, t) G is the velocity vector field of 
the gas, $(x, t) G M is the potential function of the self-gravitational force, and the stress 
^ i tensor S is given by 

5 = P/3x3-£ (^Vu + Vu*-^(divu)/3x3) -5(divu)/3x3, (1-2) 

where P is the pressure of the gas, £ > is the shear viscosity, 5 > is the bulk viscosity, 
and Vu* denotes the transpose of Vu. We consider polytropic gases for which the equation 
of state is given by 

P = P(g) = Kq\ (1.3) 

where K is an entropy constant and 7 > 1 is an adiabatic exponent. Values of 7 have their 
own physical significance pj; for instance, 7 = 5/3 corresponds to a monatomic gas, 7 = 7/5 
a diatomic gas, and 7 — )■ 1+ for heavier molecules. 
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In the simplest setting, which we consider, solutions to (11. ip are spherically symmetric. 
For r = |x|, this allows us to write 

u(x,t) = M(r,t)^ for u : [0, oo) x [0, oo) M (1.4) 

and 

gi^,t) = g{r,t). (1.5) 
The equations (11.11) then reduce to the pair 

dtQ + udrQ + ■^drijr'^u) = (1.6) 

and 

g{dtu + udru) + drP = Q{s,t)s^ds + 9, (^^^I^±^d,{r\)^ . (1.7) 

Stationary solutions q = go{r) and u = 0, which correspond to non-moving gaseous spheres 
in hydrostatic equilibrium, satisfy the following equation for Pq = Kg^: 

drPo + ^^ [ gois,t)s^ds = 0. (1.8) 
Jo 

This equation can be solved by transforming it into the well-known Lane-Emden equation 
[IJ. The solutions to (II. 8p can be characterized by the values of 7 in the following fashion 
[To] : for given finite total mass M > 0, if 7 G (6/5,2), there exists at least one compactly 
supported solution qq. For 7 G (4/3,2), every solution is compactly supported and unique. 
If 7 = 6/5, the unique solution admits an analytic expression, and it has infinite support. 
On the other hand, for 7 G (1, 6/5), there are no solutions with finite total mass. 

The stability of the Lane-Emden steady star configurations has been a question of great 
interest, and it has been conjectured by astrophysicists that stationary solutions for 7 < 4/3 
are unstable. The linear stability theory of the above stationary solutions was studied in [lOj 
in the inviscid case, namely the Euler-Poisson system, by studying the eigenvalue problem 
associated to the linearized Euler-Poisson system: any stationary solution is linearly stable 
when 7 G (4/3,2) and unstable when 7 G (1,4/3). In accordance with the linear stability 
theory, a nonlinear stability for 7 > 4/3 was established in [13] by using a variational 
approach. In the case 7 = 4/3, the analysis of [2] identified an instability in which any small 
perturbation can cause part of the system to go off to infinity. In [6], a nonlinear instability 
of the Lane-Emden steady star for 7 = 6/5 was proved based on the bootstrap argument, as 
pioneered in [4] . It is worthwhile to mention that the stability question for the Euler-Poisson 
system with 6/5 < 7 < 4/3 remains an outstanding open problem. 

The same stability question can also be asked in the presence of viscosity. There have 
been interesting studies on the stabilization effect of viscosity in the Navier-Stokes-Poisson 
system for 7 > 4/3 under various assumptions [H [H]. On the other hand, up to our 
knowledge, no rigorous stability theories are available for 7 < 4/3, the instability regime in 
the inviscid case. In this regime for viscous gaseous stars, a particularly interesting problem 
is to investigate whether or not the viscosity would dominate the gravitational force and 
stabilize the whole system. The purpose of this article is to establish the instability theory 
of the Lane-Emden steady stars whose dynamics are governed by the Navier-Stokes-Poisson 
system for 6/5 < 7 < 4/3. 

We now formulate the problem. We begin by introducing a vacuum free boundary. 
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1.1 Vacuum free boundary 

When 7 > 6/5, letting i? > be the radius of the steady star, it is well-known that 

^o(r) ~ (i? - r)^/(^-^) for r near R. (1.9) 

This boundary behavior near vacuum causes a degeneracy in the equations (11. 6p and (II. 7p . 
and it is not trivial to deal with such a degeneracy even for the local-in-time existence 
question; we refer, for instance, to [71 [TTl [12] and also [H |9] for the compressible Euler case. 
It turns out that in order to capture boundary behavior such as (II. 9p in the dynamical 
setting, one has to consider a free boundary problem associated to (II. 6p and (ll.7p as in 
[3 im [12] • We are interested in the evolution of compactly supported stars with a free 
boundary where the star meets vacuum. This is implemented by assuming there is a radius 
R = R{t) > so that 

g{r, t) > for r G [0, R{t)) and ^(i?(t), t) = 0. (1.10) 
At the free boundary we impose the kinematic condition 

j^R{t)=u{R{t),t), (1.11) 

as well as the continuity of the normal stress, Su = at the surface r = R{t). The latter 
condition reduces to 

P- — idrU--j - 6idrU + y] = for r = i?(t),t > 0. (1.12) 

Note that P{R(t),t) = Kg'^{R(t),t) = 0, so this can be reduced to a relationship between 
drU and m at r = R(t). Finally, in order for u = u{r,t):K/r to be continuous, we require 
M(0,t) = for t > 0. 

Since the boundary R{t) is free to move in time in Eulerian coordinates, it is conve- 
nient to introduce Lagrangian coordinates so that the boundary becomes fixed. Following 
the framework used in [TJ [TH [T2], we study our instability problem in Lagrangian mass 
coordinates. 

1.2 Formulation in Lagrangian mass coordinates 

We now reformulate the problem in Lagrangian mass coordinates. We set 

x{r,t) = [ A'Ks^Q{s,t)ds = [ g{y,t)dy (1.13) 

Jo J B{Q,r) 

for the mass contained in an Eulerian ball of radius r at time t. Note that 

t) = 47ir^g{r,t) (1.14) 
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and that 



dtx{r,t) = / dtg{y,t)dy = - dw{gu)dy 

'B{0,r) JB(0,r) 



/ gu ■ u = —AiTr'^g(r,t)u{r,t). (1.15) 

JdB(0,r) 



ldB{0,r) 

In particular, this imphes that dtx{R{t),t) = 0, which means that the total mass M > is 
preserved in time. The domain of x is then [0, M]. Switching to Lagrangian mass coordinates 
(x, t) G [0, M] X [0, oo) and letting the unknowns be 

p{x,t) = g{r,t) and v{x,t) = u{r,t), (1-16) 

we get the equations 

dtp + 4:7rp'^d^{r^v) = (1.17) 

and 

dtv + ATir^d^P = 167rV29^((4e/3 + S)pd^{rh)). (1.18) 
In Lagrangian coordinates our boundary conditions reduce to 

v{0,t) = 0, p(M,t) = 0, (1.19) 

and 

4:6 / V\ f 2f\ 

P - y [A'Kr'^pd^v --) -5 U'Kr'^pd^v + — \ = at x = M for alH > 0. (1.20) 
In each of these equations, we have written 

1/3 



r{x,t)=(^ r ■ (1-21) 



A simple computation, employing fll.l7p . shows that dtr{x,t) = v{x,t). 

A stationary solution p = po{x), v = 0, Pq = Kp^ to (11.171) and (11.181) satisfies the 
equation 

47irl{x)d,Po{x) + = (1.22) 

where 

dy ^ 



This is the Lagrangian version of (II. 8p . We denote such a Lane-Emden solution in Lagrangian 
mass coordinates by po with pressure Pq = Kp^. In Lagrangian x coordinates, the boundary 
behavior fll.9p is expressed as follows 

Po(a;) ~ (M - x)^/'^ for a; near M, (1.24) 

which can be also seen from fll.22p . 
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The existence and uniqueness of strong solutions to the vacuum free boundary problem 
of the Navier-Stokes-Poisson system f ll.lTp and fll.lSp featuring the behavior fll.24p of Lane- 
Emden solutions was established in [7j when 5 = 2e/3 > 0. The same methodology can be 
applied to our current setting as long as e > and 6 > 0, and we will take those strong 
solutions for granted in proving our nonlinear instability result. We remark that a well- 
posedness result in our energy space can be also proved based on our new a priori energy 
estimates for the fully nonlinear Navier-Stokes-Poisson system, described in Section HI 

1.3 Main results 

Throughout the paper we assume that 

e > 0, (5 > 0, is: > 0, and 6/5 < 7 < 4/3 (1.25) 

are all fixed. Note that although the only physical requirement on the bulk viscosity is 5 > 0, 
the assumption 5 > is critical for both our linear and nonlinear analysis. We will also write 
M, R > for the mass and radius of a stationary solution to (11.221) . 

To state the main results, we first write the system in a perturbation form. For small 
perturbed solutions a := p — po and v around the steady states satisfying (ll.22p . the Navier- 
Stokes-Poisson system (I1.17P and (ll.lSp can be written as follows: 

dta + 4:7ip'^d^{r^v) = 

dtv + Anr'd^P - Anrld^Po + ^ - 4 = lQ^'r'd,i{^e/3 + 6)pd,{rh)) ^^'^^^ 

with boundary conditions (I1.19P and (ll.20p . 

Our first main result concerns the existence of the largest growing mode of the linearized 
Navier-Stokes-Poisson system around Lane-Emden solutions, which shows a linear instability 
in the sense of Lin's stability criteria (TU] . 

Theorem 1.1. Suppose (ll.25p . There exist X > and a{x),v{x) so that a{x)e^^ and v{x)e^* 
solve the linearized Navier-Stokes-Poisson system (12. ip and (12. 2p with the linearized boundary 
conditions (12. 3 p and (12. 4p . Moreover, this growing mode yields the largest possible growth 
rate to the linearized system. 

The precise statement of Theorem 11.11 with the estimates is given in Theorem 12.11 and 
Theorem 13.21 Our second main result establishes the fully nonlinear dynamical instability 
of the Lane-Emden solutions to the Navier-Stokes-Poisson system. In the statement of the 
theorem, for any given l > and 6 > l, we write 

T':=iln-, (1.27) 
A L 

where A is the sharp linear growth rate obtained in Theorem 11.11 

Theorem 1.2. Suppose (ll.25p . There exist function spaces X and Y as well as constants 
6 > and C > such that for any sufficiently small l > 0, there exist solutions {a''{t),v''{t)) 
to (061) for t e [0, T) with T > T' so that 

\\ia^iO),v\0))\\y<CL, but sup \\ia^it),v'm\^>e. 

0<t<T' 
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The precise statement of Theorem 11.21 is given in Theorem 15.41 and the spaces X and Y 
will be clarified in Section H] and |3 



Remark 1.3. Our results show that regardless of how large the viscosity parameters e, 6 are, 
and no matter how small smooth initial perturbed data are taken to be, there is no stabilization 
of the system. We conclude from this that all Lane-Emden steady star configurations for 
6/5 < 7 < 4/3 are unstable, regardless of viscosity. 

Remark 1.4. The escape time T' is determined through (11.271) by the linear growth rate A. 
We note that the instability occurs before the possible breakdown or any collapse of strong so- 
lutions. We also remark that the instability occurs in the X norm, which can be characterized 
by the instantaneous physical energy 



M / ^ 
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-f)|up H ~P] [ ^p^^^^c/x(iy in Eulerian coordinates. 

Of course, this is not a coincidence: the Lane-Emden solutions for 7 < 4/3 do not minimize 
the physical energy functional and thus one might expect some kind of instability. Note that 
they do minimize for 7 > 4/3 (see for instance JBi H^). 

The presence of viscosity and the nonlinear boundary condition f ll.20p for the Navier- 
Stokes-Poisson system make the problem distinguishable and interesting not only from a 
physical point of view, but also from a mathematical point of view. What follows now are 
some of the main mathematical difficulties we encounter in analyzing the system, and a brief 
discussion of our methods for resolving them. 

The proof of Theorem 11.11 is based on a variational analysis of equations obtained by 
linearizing fll.26p . The main difficulty that arises in constructing growing- mode solutions is 
that, due to the viscous terms, the growth rate (eigenvalue) appears in the problem with 
two different homogeneities. This breaks the natural variational structure used in [10] to 
construct growing modes in the inviscid case. To get around this difficulty, we employ a 
technique introduced in [5]: we introduce a relaxed parameter that allows us to remove 
one of the eigenvalue homogeneities, study the resulting modified eigenvalue problem (which 
has a nice variational structure), and finally return to the original formulation through a 
fixed point argument. While the solutions constructed in this manner are definitely growing 
modes, it is not clear a priori that they grow at the largest possible rate. To verify this, 
we carry out a careful analysis, paying particular attention to the boundary behavior of the 
growing mode, which will be crucially used in the subsequent nonlinear bootstrap argument. 

The proof of Theorem 11.21 is based on a bootstrap argument from linear instability to 
nonlinear dynamical instability. Passing from a linearized instability to nonlinear instability 
requires much effort in the PDE context since the spectrum of the linear part is fairly 
complicated and the unboundedness of the nonlinear part usually yields a loss in derivatives. 
In order to get around these difficulties and to find the right space Y, we employ careful 
nonlinear energy estimates for the whole system so that: first, the nonlinear estimates can be 
closed; and second, their interplay with the linear analysis can complete the argument. For 
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this particular problem, the space Y is minimally chosen so that the viscosity disturbance 
near the vacuum boundary can be controlled within Y. We note that in Lagrangian mass 
coordinates, the continuity equation interacts well with the viscosity term, which allows us 
to derive nice estimates for a/po and its temporal and spatial derivatives. This plays an 
important role in closing our nonlinear energy estimates. 

Another delicate and important issue is the nonlinear boundary condition f ll.20p . In 
order to carry out higher-order energy estimates that require integration by parts, we can 
only employ differential operators that respect the boundary conditions, namely temporal 
derivatives. This forces us to carefully use the structure of the equations in order to gain 
bounds on spatial derivatives. A second difficulty with the boundary arises because we use 
Duhamel's principle to study the nonlinear problem with the linearized evolution operator. 
The linearized boundary condition is homogeneous, but the nonlinear boundary condition 
is certainly not. This forces us to introduce a corrector function that removes the boundary 
inhomogeneity. While the construction of this function is not particularly delicate, the 
regularity required to do so dictates that we close our energy estimates at a higher order 
than we would otherwise. 

The paper proceeds as follows. The first half is devoted to the development of the linear 
theory and the proof of Theorem 11.11 In Section we formulate a variational problem to 
find a growing-mode solution to the linearized Navier-Stokes-Poisson system. In Section [3l 
we show that our growing-mode solution grows at the largest possible rate. In the second 
half of the paper, we carry out our nonlinear analysis. In Section 4 we derive high-order 
nonlinear energy inequalities. Based on the linear growth and the nonlinear estimates, we 
then prove the bootstrap argument and Theorem 11.21 in Section 5. 

2 Construction of a growing mode solution to the lin- 
earized equations 

2.1 Linearization around a stationary solution 

We now linearize the equations in Lagrangian mass coordinates around the stationary solu- 
tion V = 0, p = po, r = ro (as defined by f ll.23p ). We will write a for the linearized density, 
and (by abuse of notation) v for the linearized velocity. Then the linearized equations are 
given by 





and 



Ti-rg Jo Po{y) 

where we have written P = '-fKpQ~^a. The linearized boundary conditions are 

v{0,t) = 0, a{M,t) = 0, 




(2.2) 



(2.3) 



and 



P - — Anrlpod^v 




6 inrlpod^v + 




at X = M for alH > 



(2.4) 
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Again, we can view (12. 4p as a boundary condition only for v since P = jKpl'^a = at 
X = M for each t > 0. 

It will often be useful for us to analyze a variant of this system, where we analyze the 
unknowns a and w := TqV. For these unknowns, the equations fl2.1l) - fl2.4p become 

dta + Anpldxiv = 

dtw + Airr^d^i'yKp^'^a) - Avod^Po ^^^dy = levrV^a^, + 6^ pod^w ^ 

along with the boundary conditions 

-2 (0, t) = a{M, t) = and — Airrlpod^ [ — ]] + SiAnpod^^w) = at x = M. (2.6) 



2.2 Growing mode solution 

We want to construct a growing mode solution to the linearized equations. We do so by 
looking for a solution of the form 

a(x, t) = cT(x)e^* and v{x, t) = v{x)e^^ (2.7) 

for some A > 0. If we can find such a solution, then we say the solution is a growing mode 
since |e'*'*| — )• oo as t — t- oo. Plugging the ansatz (12. 7p into the linearized equations (l2.ip -f lT^ 
and eliminating the time exponentials, we arrive at a pair of equations for (t{x) and v{x): 

Xa + ATipldxirlv) = (2.8) 

and 

Xv + Anrld^P + ^ [ 4r\'^y = 167rVo'9,((4£/3 + 6) pod,{rlv)) , (2.9) 
along with boundary conditions 

v{0) = a{M) = and - y (^4:7rrlpod:,v - - 5 (^ATcrlpod^v + — ^ = at a; = M. 

(2.10) 

Our main result of this section establishes the existence of such a growing mode. 

Theorem 2.1. There exist X > and a,v : (0, M) R that solve f l2:8|) - (l230|) and satisfy 
the following. 

1. a and V are smooth on (0, M) and satisfy (I2.8p - (l2.9p classically for x G (0,M). 

2. It holds that 

\'v{oc]\ 

limsup — -— - + limsup |(t(x)| + limsup (9i:(rQt')(x) < oo. (2-11) 

a;-s>0 rQ[x) x-^0 x-s>0 

In particular, t>(0) = 0. 
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3. LetD denote the linear operator D = pod^- ThenTi^v andD^{a/pQ) have well-defined 
traces at x = M for every integer k > 0. In particular, a{M) = 0. 

4- A > satisfies the following variational characterization: 

2^ 



A 



4e 



-Po 



40^ 



9 







dx 



+ 



M 



V 2 ' ' 2Tirl 



^ ' 'e\' ] dx > -X' 



M 



dx (2.12) 



for every 6 satisfying y/p^d,J G L^{{0, M)) and 6/{rl^) G ^^((O, M)). Note that for 
such 6, it holds that 6'/(rQ^/po) G L^((0,M)), which means that all of the integrals in 
fl2.12p are well-defined. 

5. It holds that 



M 



Po 



Po 



dx 




+ Po \dx{rQv)\ + rg \dx{podx{rQv))\ ] dx < oo. (2.13) 



The proof of Theorem 12.11 will be completed in Section 12.61 Throughout the rest of 
the section we develop the tools needed in the proof. First we reformulate fl2.8p -( |2rT0l) 
to involve a single unknown function, (p. The resulting problem for (p does not possess a 
standard variational structure since A appears both linearly and quadratically. To construct 
a solution using variational methods (required for proving fl2.12p . which is essential for the 
linear estimates of Section |3]), we employ the technique of Guo-Tice [5], which proceeds as 
follows. We modify the problem by replacing the linear appearance of A by an arbitrary 
parameter s > 0. The resulting family (every s > 0) of problems is amendable to solution 
by the constrained minimization of an energy functional, and for a range of s we show that 
A = A(s) > 0. We then study the behavior of A(s) as a function of s and show that it is 
possible to find a unique fixed point so that A(s) = s > 0. This then yields the desired 
solution (p, which in turn yields the solution to ( l2.8p -( l2TT0l) . 

We begin by reducing to the study of a single unknown by introducing the function 



We may then use fl2.8p - fl2.10p to compute 

A 



47rrQ 



cr = Po^^, and d^P = d^i^ypoPod^ 



(2.14) 



(2.15) 



where Pq = Kp^. Using these and replacing in (12. 9p . we arrive at a second-order equation 
for 6: 



d. 



4Xe 



3 / / nrn 



+ X5 + 7P0 Podx'P + 



dxPo 



X' 



IGvr^rn 



(2.16) 
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The corresponding boundary conditions are 

^(0) = and yA (^47rr^po<9x + ^A(47rpo<9x0) = at x = M. (2.17) 

2.3 Modification of the problem 

Note that Theorem 12.11 is phrased in Lagrangian mass coordinates. This is because we will 
use these coordinates in our nonlinear analysis later in the paper. However, constructing 
the solution to f l2.16p -f l2.17p is somewhat easier if we make a change of variables back to 
the Eulerian radial coordinates associated to the stationary solution. To avoid confusion 
with the Eulerian radial coordinate for the nonlinear problem, we will call our new variable 
z = ro(x), where vq is given by f ll.23p . If a; G (0, M) for M the mass of the stationary star, 
then z G (0, R) for i? > its radius. We will write go{z) = po{x) for the stationary density, 
Pq = Kg'^, and (p{z) = (f){x) for the new unknown in z coordinates. Then 

= -r^d,. (2.18) 

In these coordinates, the equations fl2.16p becomes 

For the boundary condition at 2; = i? we use fl2.17p to see that 

M^M?),1^(^?^_3^).0. (2.0) 

At 2; = we enforce the boundary condition V5(0) = 0. Once we have a solution in hand, 
we will show that, in fact, ip{z)/z'^ — as 2; — )■ 0, which allows us to switch back to the 
boundary condition (0/rg)(O) = 0. 

There is a difficulty in viewing f l2.19p -f l2.20p in a variational or Sturm-Liouville framework 
because of the appearance of A with two different homogeneities. To get around this issue, 
we temporarily modify the problem in order to restore the variational structure. Ultimately 
we will undo the modification and return to the proper formulation. 

Fix s > and define 

e = se and 5 = s5. (2.21) 
Instead of 02.191] . we will analyze the equation 

- 5. ^ (^y + 5 + 7P0 j — J + 4y^^ = --^^ (2.22) 
for arbitrary s > 0. We couple this equation to the boundary conditions V'(O) = and 

Modifying the problem in this way restores the variational structure. Indeed, in fl2.22p 
the A^ term can be viewed as an eigenvalue. Thinking of the principal eigenvalue A as a 
function of s, i.e. A = A(s), we will show that it is possible to choose s so that A(s) > and 
s = A(s), which returns us to the original problem and yields a growing- mode solution. 
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2.4 Constrained minimization 

In order to construct solutions to f l2.22l) - fl2.23l) . we will employ a constrained minimization. 
To begin, we define the function space on which the energy functionals will be defined. 
For r > we define the weighted Sobolev space H^{{0,R)) as the completion of {u G 
C°°{[0,R]) I u{0) = 0} with respect to the norm 



Ml^ = I ; dz, (2.24) 

Jo ^ 

where ' = d/dz. This weighted Sobolev space possesses the same sort of embedding (con- 
tinuous and compact) properties as the usual space H^. Since these results embeddings are 
not widely available in the literature, we record them in the following lemma. 

Lemma 2.2. The following hold. 

1. For u G H^{{0, R)) , we have the inequalities 

sup k(z)z-(^+i)/2| < ( Mfl^rf^^ (2.25) 

0<z<R 

and 



u'iz)\^ + \u{z)\ 



«K^)|^_ 4 f^\u'iz)\ 



,-dz<- -/ ' ' ^' dz. (2.26) 

2. Let < a < 1. We have the compact embedding Hl{{0,R)) CC L^_,_]^_,_„((0, i?)), where 
the latter space is the weighted space with norm 

Ml. = r^-^^dz. (2.27) 

Proof. We begin with the inequalities in item 1. By approximation we may assume that u 
is smooth and u{0) = 0. Then 

1/2 / ,^/(^)|2 N^/^ 



\u{z)\ = \u{z) - u{0)\ < I \u'{t)\dt<{ I t^dt] I / ' dt 

^(S)"(f^*)"' 

which yields the first inequality. To get the second, we recall an inequality due to G. H. 
Hardy: 

for 1 < p < oo and < 6 < oo. Then \u{z)\ < \u'{t) \ dt implies that 

^ |m(2)P , f f\ n .V dz 

' ' dz< / \u'{t)\dt) --. (2.30) 



^ Jo \Jo J ^ 
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Applying Hardy's inequality to the right side with / = u'x{o,r), b = t + 1, and p = 2 yields 

i"Wi!,,<^rM££,,, ,2.31) 



z-+^ - (r + 1)2 Jo 

which is the desired inequality. 

We now prove the compactness result. Assume that ||^in|li^i < C for G N. Fix k > 0. 
We claim that there exists a subsequence {ur^} so that 

sup - M„H 2 < K. (2.32) 
To prove the claim, let Zq & {0, R) be chosen so that 

Then since the subinterval {zo,R) avoids the singularity of l/z'^, Un\{zo,R) is uniformly 
bounded in H^{{zo, R)). By the compact embedding H^{(zo, R)) CC C^{{zo, R)) we may 
extract a subsequence {wn,} that converges in L°°{{zo, R)). We are free to restrict the sub- 
sequence to sufficiently large values of i so that 

Then along this subsequence we can apply the first inequality in item 1 to get 

^ I Un, (z) - Un. {z)? , /""^f I (z) - M„, (2;) I \ \ Un, {z) - n„ . {z) I ^ 



which proves the claim. Now we may use the claim with k = 1/k, k & N and employ a 
standard diagonal argument to extract a subsequence converging in L'^^^^^{{0, R)). □ 

Remark 2.3. The inequality fl2.26p implies that we can take the norm on to be 

llu\\l,= n.!^^dz. (2.36) 



We can now define the energy functionals to use in the constrained minimization. Let 



R , .in _|2 



^2 Sz'^ 



z 



2 .aPf 



E{^)= I (5 + 7Po)^ + Tf^ 5.^-3^ +4^l¥'l' (2.37) 



z^ 



and 

^(¥')= r^lv'l'- (2-38) 
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Z2 



By (I2.26P in Lemma |2.2[ both E and J are well-defined on the space Hl{{0, R)). Note, 
though, that E is not positive definite since d^Po < 0. Define the set 



A:= {^eH',{{0,R))\J{^) = l}. 



(2.39) 



We will build solutions to f l2.22p by minimizing E over A. First we show that such a 
minimizer exists. 

Proposition 2.4. E achieves its infimum on the set A. 

Proof. To begin, we show that E is coercive on A, which amounts to controlling the last 
term in E. Recall that by (11.91) . go{z) ~ (i? — zY^^"'~^'^ for z near R. This implies that 



Qo 



-fKgl d;,go = rdziQa ) 

7-1 



(2.40) 



is bounded near z = R. Since Qq and Pq = are smooth and bounded below away from 
z = R, this implies that 

d.Po 



Qo 



Then for any zq E {0, R) we may bound 



< oo. 



(2.41) 



L°°{{0,R)) 



R 



\d.Po 



m 
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\d.Po 



< Zo II -Poll L° 



+ 



R 



\dzPo\ Qo Iv 
ZQo 



r 'rf - 


1 

h — 


dzPo 


QoW? 


Jo z 


Zo 


Qo 


1 



< Zo-^ \dzPo\\ L° 



R 



\dz'^\ 



1 

+ — 

Zo 



dzPo 



Qo 



(2.42) 



For the second inequality we have used Lemma 12.21 and the fact that ip & A. Then by 
choosing zq sufficiently small, we have that 



R 



6 



+ 7^0 



+ 



3^ 



dz^ - 3 



(2.43) 



for a constant Cz > depending on the choice of z, which immediately yields the desired 
coercivity since 6 > 0. 

With the coercivity in hand, we may deduce the existence of a minimizer by using the 
standard direct methods, employing Lemma 12.21 for compactness. 



Since a minimizer exists, we can now define the function fi : (0, oo) — > M by 

/i(s) = inf E{(p; s), 



□ 



(2.44) 
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where we have written E{ip) = E{ip; s) to emphasize the dependence of E on the parameter 
s > 0, i.e. 



E{^;s) 



R 



\dz<f\ 



+ 



Ae 
3? 



+ / iPo 



\dzV\ 



+ 4- 



(2.45) 



The minimizer we have constructed satisfies Euler-Lagrange equations of the form f l2.22p . 



Proposition 2.5. Let E A be the minimizer of E constructed in Proposition 2.J^. Let 
fi := E{(f). Then (f is smooth on {0,R] and satisfies 



dz^p 



(2.46) 



along with the boundary conditions ip{0) = and 
Proof. Fix ipo eH^i{0,R)). Define 

j(t,T) = J{(p + tipo + Tip) 

and note that j(0, 0) = 1. Moreover, j is smooth and 



|(0.0) = 2 



R 



dj 



and ^(0,0) = 2 / go^ = 2 



(2.47) 



(2.48) 



(2.49) 



So, by the inverse function theorem, we can solve for r = r(t) in a neighborhood of as a 
function of t so that r(0) = and j{t, r(t)) = 1. We may differentiate the last equation 
to find 

%(0,0) + |^(0,0)r'(0) = 0. 



dt 



Or 



and hence that 



r'(0)^-i|(0.0)^-^ 
Since if is a. minimizer over A, we then have 



R 



Qo- 







di 



E{ip + tipQ + T{t)ip), 



(2.50) 
(2.51) 

(2.52) 



which implies that 

'''~5 + lP^ 











(/.(y^o + r'(0)</.) 



+ 



^ 4£ 
3^2 
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Rearranging and plugging in the value of t'(0), we may rewrite this equation as 




where the eigenvalue is /x = E{ip). 

By making variations with cpo compactly supported in (0, R) we find that if satisfies the 
equation fl2.46p in a weak sense in {0,R). Standard bootstrapping arguments then show 
that ip G H^{{zq, R)) for all /c > and < < -R, and hence if is smooth in (0,-R]. This 
implies that the equations are also classically satisfied. Since tp G H'^{{R/2, R)), the traces of 
(p, dzif are well-defined at the endpoint z = R. Making variations with respect to arbitrary 
V?o £ C^{{0, R]), we find that the boundary condition f l2.47p is satisfied. The condition 
(p{0) = is satisfied by virtue of Lemma [2^ 

□ 

We now want to show that the minimizers, which are solutions to fl2.46p . satisfy the 
asymptotic condition |v5(-z)| /-z^ — as 2; — )■ 0. As a preliminary step we record an asymptotic 
result for solutions to a more generic ODE. 

Lemma 2.6 (Proposition A.l of [lO]). Suppose that iP{t) solves 

i)"{T) + {ar-^ + ^(r))^'(r) + r-V(r)^(r) = on < r < Tq (2.55) 
where ' = d/dr, and f,g G C°([0,ro]). If a < 0, then either iIj{0) ^ or 

mr)\ < and \^'{r)\ < ^. (2.56) 

Proof. The case a > 2 is the content of Proposition A.l of [10], but the proof of the propo- 
sition also shows the result when a < 0. □ 

Next we use this lemma to establish the asymptotics at 2; = for solutions to f l2.46p . 

Lemma 2.7. Suppose ip is a solution of fl2.46p . Then \'p{z)\ < Cz^ and \ip'{z)\ < Cz^ near 
2 = 0. 

Proof. We begin by rewriting fl2.46p . Define X = e + 5 + 7P0 = e + S + 'yKg^ and Xq = 
e + Pq = e + KqI. Then fl2.46p is equivalent to the equation 

Note that X'/X, X^/X, and Qq/X are all continuous at 2; = 0, so we may apply Lemma 
12.61 with a = — 2 to deduce that either ip{0) 7^ or \(p{z)\ < Cz^ and |v2'(-2^)l ^ Cz"^ near 
z = Q. By Lemma 12.21 the former condition cannot hold, so the latter conditions must be 
the case. 

□ 
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2.5 Properties of the eigenvalue //(s) 

We are ultimately concerned with finding /i — — for some A > 0. This requires us to work 
in a range of s so that ii{s) < 0. Our next result shows that //(s) < for s sufficiently small. 

Lemma 2.8. There exist constants Ci,C2 > depending on 7, so that 

Ks) < sCi - C2. (2.58) 

In particular, < for s sufficiently small. 

Proof. Via the homogeneity of E and J we may write 

l,{s) = M{E{ip; s)/J{ip) I cp e H^iiO, R))}. (2.59) 

The structure of E points to a natural choice for a test function: (p{z) = z^. Using this, 
integration by parts reveals that 

E{ip;s) = J^ (s + -fPo^9z^ + Ad,Poz^ = 95 s z'' + - 12) Poz\ (2.60) 

Since 7 < 4/3 it holds that (97 - 12) < 0. This implies that 

Ks) < = sC, - C2, (2.61) 

where ^ 

Ci = {3SR^) (^j^ go{z)z^dz^ > (2.62) 

C2= (^J^ {12- 9-f)Kgl{z)z^dz^ (^J^ go{z)z^dz^ > 0. (2.63) 

□ 

The next proposition proves some crucial monotonicity and continuity properties of ii{s) 
for s > 0. 

Proposition 2.9. The following hold. 

1. ii{s) is strictly increasing in s. 

2. There exist constants C3, C4 > 50 that 

H{s) > -C3 + SC4. (2.64) 

3. n E Ci^{{0, 00)), and in particular /j, e C°((0, 00)). 



and 
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Proof. We begin by establishing two bits of notation. According to Proposition \2A\ for each 
s G (0, oo) we can find (/^s G ^ so that 

E{ifa; s) = inf E{(f; s) = fi{s). (2.65) 

Next, we decompose E according to 

E{^; s) = Eo{^) + sEiiy^) (2.66) 

for 

Eo{^):= [%Po^^ + ^^\^f (2.67) 
Jo z z 

and 

EAv)--= r^^ + ^ ^^^-3-'>0. (2.68) 



The non-negativity of Ei imphes that E is non-decreasing in s with ip & A kept fixed. 

To prove the first assertion, note that if Si, S2 G (0, 00) with Si < Si then the minimahty 
of ipsi and the non-negativity of Ei imply that 

/i(si) = E{ip,^\ si) < E{ips^\ si) < E{lPs2; S2) = yu(s2)- (2.69) 

This shows that /i is non-decreasing in s. Suppose by way of contradiction that /i(si) = /i(s2). 
Then the last inequality implies that 

siEiips,) = S2Ei{ips,), (2.70) 

which means that Ei^ips.^) = 0. The vanishing of Ei^ips^) implies that ips2 = 0? which is 
impossible since (ps2 ^ Hence equality cannot be achieved, and /i is strictly increasing in 
s. 

Now note that ( 12. 66^ and the non-negativity of Ei imply that 

/i(s) > inf Eo((^) + s inf Ei((^). (2.71) 



The proof of Lemma 12.81 shows that 



and it is a simple matter to see 



:= inf Eo{p) < 0, (2.72) 

ip&A 



C4 := inf Ei(v9) > 0. (2.73) 

•f&A 



The second assertion follows. 

Now fix Q = [a, 6] CC (0, 00), and fix any ip E A. Again by the non-negativity of Ei and 
the minimality of ips we deduce that 

£;(V'; b) > E{ilj- s) > Eips; s) > aEi{ps) - C3 (2.74) 
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for all s & Q. This implies that there exists a constant < C = C{a,b,ip, 7, K) < 00 so that 

supEi{ifis)<C. (2.75) 

Let si, S2 G Q. Using the minimality of ips^ compared to ips2, we know that 

= E{(ps^; si) < E{ips2; si), (2.76) 

but from our decomposition fl2.66p . we may bound 

E(ips2; si) < E{(fs2\ S2) + \si - S2I Ei{ips2) = /i(s2) + |si - S2I Ei^ips^). (2.77) 

Chaining these two inequalities together and employing fl2.75p . we find that 

/i(si) </i(s2) + C|si-S2|. (2.78) 

Reversing the role of the indices 1 and 2 in the derivation of this inequality gives the same 
bound with si switched with S2- We deduce that 

- Ks2)\ < C \si ~ S2\ , (2.79) 

which proves item 3. 

□ 

Now we know that the eigenvalue fi{s) is negative as long as s < C2/C1 and that fi is 
continuous on (0, 00). We can then define the non-empty open set 

= /i-^((-oo,0)) C (0,cx)), (2.80) 

on which we can calculate A(s) = ^y—fi{s) > 0. 

It turns out that the set Q is sufficiently large to find s > so that A(s) = s. This 
inversion will then allow us to solve the original growing-mode equations. 

Proposition 2.10. There exists a unique s E so that X(s) = a/— /i(s) > and X{s) = s. 

Proof. According to Propositions 12.81 and 12. 9| we know that fi{s) < for s G [0,C2/Ci). 
Moreover, the lower bound fl2.64p implies that fi{s) — )■ +00 as s — ?■ 00. This implies the 
existence of sq G (0, 00) so that Q = (0, sq), which means that A(so) = 0. Define the function 
\Ef : (0, So) — (0, 00) by \E'(s) = s/A(s). The monotonicity and continuity properties of /i are 
inherited by \E', i.e. \E' is continuous on (0, sq) and strictly increases from to +00 as s — )■ sq. 
As such, we may apply the intermediate value theorem to find a unique s G (0, sq) so that 
\E'(s) = 1. For this s we then have that s = A(s), the desired result. □ 
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2.6 Proof of Theorem [231 

We now combine our above analysis to deduce the existence of a solution ip, X > to 

Theorem 2.11. There exists A > and (f G H2{{0, R)), smooth on (0, R], that solves f l2.19p 
along with the boundary condition (I2.20p . The solution satisfies the asymptotics \^{z)\ < Cz^ 
and \dz(f{z)\ < Cz^ as 2; — )• 0. 

Proof. Combining Proposition 12.51 and Proposition I2.10[ we see that there exists a solution 
to fl2.22p and fl2.23p for A(s) = a/— > 0, satisfying s = A(s). This imphes that the 
solution is actually a solution to fl2.19p and ( I2.20p . The asymptotics at z = follow from 
Lemma 12.71 □ 

Remark 2.12. It is actually possible to get a lower bound for the size of X > in terms of 
6, Qo, 7, and R. Indeed, from Lemma \2. <^ it holds that X^ + ACi — C2 > 0, and hence 



A > -^' + f™ > 0. (2,81) 

Here the constants Ci,C2 > are determined explicitly in terms of 6, go,'~f, and R according 
to fl232|) and fl233D . 

An immediate consequence of Theorem 12. Ill is the existence of a solution to fl2.16p - fl2.17p . 

Corollary 2.13. There exists A > and (f){x) = ip{ro{x)), smooth on (0,M), that solves 
fi2.l6p -f l2.i7p . The solution satisfies 

\6(x)\ 

limsup ^VHr + hmsup \d.Mx)\ < 00. (2.82) 
x^o r^{x) x^o 

Let T) denote the linear operator D(f){x) = po{x)dx(f){x) . The solution satisfies the property 
that D^cj) has a well-defined trace at x = M for every integer k > 0. 

Proof. All of the conclusions, except those concerning D follow directly from Theorem 12.111 
When k = 0, the trace DV(M) = 0(M) is well-defined since f{R) = <y5(ro(M)) is well- 
defined. Note that 

m^) = Po(x)9.^(x) = ^ mM) = (2.83) 

SO that 2)(/)(M) is well-defined. We may argue similarly, using the fact that d^ip{R) is well- 
defined for all > 0, to deduce that 'D^(j){M) is well-defined for all > as well. 

□ 

Now, with Corollary 12.131 in hand, we are ready to present the 
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Proof of Theorem \2.1[ Let A > and be the solution to fl2.16p -( l2.17p given in by 
Corollary I2.13[ Let us then define v and a according to 

A 



2 



and a = pld^cj). (2.84) 

' 

Using these definitions of v and a in conjunction with the properties of recorded in Corollary 
I2.13[ we easily deduce the first three items of the theorem. 

To prove the variational characterization of the fourth item, we return to the variational 
characterization of A in z = ro(x) coordinates. According to Theorem 12. IH A > satisfies 



> -A' I %W\^dz (2.85) 
Jo ^ 

for every 'd G H\{(^^I{)). Then the variational characterization in fl2.12p follows by making 
a change of coordinates Q{x) = 'dlz) = ?9(ro(x)). Note that i} G H^ilOjR)) if and only if 
G L2((0,M)) and 0/{r^^) G L^{{0,M)). Also, changing coordinates in (E^SD of 
Lemma 12.21 shows that ^/(?"oa/Po) ^ -^^((0,Af)), which means that all of the integrals in 
(12J2|) are well-defined. 

We now turn to the proof of fl2.13p . Using the inclusion ip G H2{{0, R)), the above 
analysis implies that y/f)odx4>, (p/ {r^y/po) G L^((0,M)). From this and equation f l2.16p we 
may then deduce that 



/oPo 

This and fITMD then imply fl^TT^ . 



/ \^ + Po 1^x01' + ^0 mPodM'] dx < cx). (2.86) 
^0 V '^oPo / 



□ 



3 Linear estimates 

Due to the indirect way in which we constructed growing mode solutions in Section [2|, it is 
not immediately obvious that the A > of Theorem 12.11 is the largest possible growth rate. 
However, because of the inequality ( 12.12p . we can show that no solution to the linearized 
problem f l2.ip -( l27il) can grow in time at a rate faster than e^*. Hence the growing mode 
constructed in Theorem 12.11 actually does grow in time at the fastest possible rate. The proof 
of this result and its implications for solutions to the inhomogeneous linearized problem are 
our the subject of this section. 

3.1 Estimates in the second-order formulation 

First, we will prove estimates for solutions to the following second-order problem. 

forxG(0,M) (3.1) 



-1^ 2 4 = - Ox 

lOTT^rn TCTn 



4e 

— + podxdt(l) + -fPopodx 



3 
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with boundary conditions 



0(0, t) = and — Anr^pod^ ^ 







(5(47rpo40) = at X = M, 



(3.2) 



and initial conditions (j){x, 0) and dt<f){x, 0) given. We will assume throughout that (p satisfies 
Vp^9,0 e L2((0,M)) and ^/(r^yp^) G L2((0,M)). 

Solutions to this linear problem obey an energy evolution equation related to the inequal- 
ity f l2.12p . We record this now. 

Proposition 3.1. Suppose (p is a solution to fl3.ip - (l3.2p . Then 



327r2r, 



-jdx + [ \5po + ^Po r^d-^ ( ^ 

Jo \ 3 V ^0 



dx 



+ 



27rrf 



|0r dx. (3.3) 



Proof. Multiply (13. ip by ^0 and integrate over x G (0,M). An integration by parts, an 
application of the boundary conditions (13. 2p . and some simple algebra yield the desired 
equality. □ 

We can use this and the variational characterization of A given in Theorem 12.11 to deduce 
some estimates of the following norms: 



2 . 

1 •" 



M 



I0f 



IGvr^rQ 



dx 



2 

2 



i^po \dx(t)\^ + yPo rldx {^^^ j dx 



2 . 

3 



M 



Our estimates for these norms are the subject of the next result. 
Theorem 3.2. Let cf) solve (I3.ip - (l3.2p . Then we have the following estimates: 

\2 



(3.4) 
(3.5) 
(3.6) 



^ + ^ ||0(.)||^c^.<e^^* 110(0)11^ + §(e^^*-l). 



1 
A 



ia0(t)ll? + ||0(t)ll^<e2^* (2A 



and 



Here 



and 



110(0)11? + §(e^^*-l) 



la 



1 + 2 



2 27^0 



X 



Co = 2 sup 3 

xG{0,M) ^ol^J 



A 



< oo. 



2 

2 ' 



(3.7) 
(3.8) 

(3.9) 

(3.10) 
(3.11) 
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Proof. We integrate the result of Proposition 13.11 in time from to t to see that 



A/ 



327r2r^ 



dx 



t rM 



^0 

27rrQ 



dxds 



' dx — 



M 



iPopo 



\dMt)r + 7^mt)\ ]dx. (3.12) 



Note that since dxPo = — x/(47rrQ), we have that 



' dx 



X 



dx < 0. 



(3.13) 



The variational characterization of A given in fl2.12p of Theorem 12.11 allows us to estimate 



-U{ 



\dA{t)? + ^^\m\']dx 



M 



6po \d^(f){t)f + yPo 



' 



dx < — 



2 Jo lQn^4 



dx. (3.14) 



We may then combine f l3.12p -( l3.14p to see that 



A/ 



<Ko + 



dx 



t rM 



^0 



5po \dxdt(t){s)\^ + yPo 



rid. 



,2 rM 



levrV^ 2 



M 



Spo \dx4>{t)f + yPo 



rid. 



dxds 
rl J 



dx, (3.15) 



which we may rewrite as 

lmm\l + l\\dt<P{s)\\ids<Ko + ^ 

Integrating in time and using Cauchy's inequality, we may bound 



(3.16) 



A||(/)(t)||^ = A1|0(O)1|^ + A / 2{<p{s),dt<P{s)hds 



< A 110(0)11^+ / \\dt(j){s)\\'^ds + y 



s)\\lds. (3.17) 



On the other hand 



Xdt \\m\i = A2(9,0(t),0(t))i < A^ \\m\i + \\dt 
We may combine these two inequalities with (13.160 to derive the differential inequality 



(3.18) 



dt\\m\i + \\m\2<Ki+2x 

22 



t + 2X 



s)\\^ds 



(3.19) 



for Ki as defined in the hypotheses. An apphcation of Gronwall's lemma then shows that 



+ 



s) ds < e 



2Xt 







(3.20) 



for all t > 0, which is the bound (13.71) . 

To derive the estimate (13. 8p . we return to (I3.16P and plug in (I3.17P and (I3.20p to see that 



j\\dtm\\l+\\m\l<K,+x 



1 + 2A 



s)\\lds<e'''{2\\\m\\l + Ki). (3.21) 



Finally, for (13. 9p we return to (I3.12p and employ (13.131) to see that 



, ,3 < '<« - 

Since L'Hospital's theorem implies that 



X 



rim' 



X 4:iTpo{x) 47rpo(0) ^ 

lim = < 00, 



hm 



we may deduce that 



X 



(0,M) ^ol^J 



< 00. 



(3.22) 



(3.23) 



(3.24) 



The estimate (13. 9p then follows directly from (]3.22p . (]3.24p . and the estimate of ||0(t)||i in 



□ 



3.2 Estimates in the first-order formulation 



Now consider a and w to be solutions to the first-order linear system (12.51) with boundary 
conditions (12.60 and initial conditions cr(a;, 0) and w{x, 0). A simple calculation shows that if 
we apply dt to the second equation in (12. 5p and then eliminate dfCr by using the first equation 
in (12. 5p . then we arrive at the second-order formulation (I3.ip - (l3.2p for (p = w. Then Theorem 
13.21 yields various estimates for w = (j). We now seek to rewrite these estimates for w and to 
use them to derive a similar estimate for cr. 

Theorem 3.3. Let a,w solve the linear system (I2.5p - (l2.6p . Then 



M 



W{t)\' 



■dx + 



M 



,7-1 




<j{t) 



Po 



dx 



+ IQn' 1^'' i^6po \dMt)\' + f Po rid. (^^^ 



Po 



dx < Ce^^* 
\dtw{0)\ 



M 



■dx 



dx 



M 



6po \d.w 



2^45 

+ yPo 



w{0) 



dx 



(3.25) 
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Proof. We switch to the second-order formulation for (j) = w. Then the estimates (j3T])- (13. 9p 
of Theorem 13.21 imply that 



Jo ^0 Jo 



+ 167r' 



AI 



Spo \dxw{t)f + ypo 



3n / 



dx < Ce 



2Xt 



M 



■dx 



' 



+ 



M 



+ 



M 



Spo \dxW{0)f + ypo 



^O'-'x \ 3 



dx 



. (3.26) 



Let us call the term in the brackets on the right side of this equation Zq. Since dta 
—4:TTpldxW, we then have that 



M 



jKp2+'\dMt)\"dx 



IGtt'^pI 



\dt(j{t)f dx := \\dt(J 



We then have that 



dt\\ait)\\,<\\dMt)h<VCZ',e^'. 
Integrating this in time, we then find that 



kWll4< 11^(0)114 + 



A 



e^*-l) <Ce^V 11^(0)114 + ^0- 



The estimate (Km then follows directly from ([326]), (Km . and (Km . 

3.3 Estimates for the inhomogeneous first-order problem 

Consider the linear operators 



Ciw = inpld^w 
C^a = 47rr^9^(7irp2;"V) - Arod^Po 

Jo 



Poiy) 



C3W = -1677^0 (9^ 

and the corresponding matrix of operators 

C-- 

We also consider the boundary operator 



4e 



Ci 
C2 Cs 



+ 6 pod^w 



B{w) = -y (^iTir^podx (^)) ~ ^(^^'^o^^ 



w). 



(3.27) 
(3.28) 

(3.29) 
□ 



(3.30) 
(3.31) 

(3.32) 
(3.33) 

(3.34) 
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Notice that the first-order equations f l2.5p -( l2l6|) are equivalent to the equation 



with homogeneous boundary conditions 



w 

— (0, t) = a{M, t) = and B{w) =0 at x = M. (3.36) 



Let us denote e*^ for the solution operator to fl3.35l) - fl3.36l) . i.e. 



where a and w solve fl3.35p - fl3.36p with initial data cr(0) and w{0). 
Suppose now that cr and w solve the inhomogeneous problem 



.2 

along with the boundary conditions 

w 

^(0, t) = a{M, t) = and B{w) = Nq at x = M. (3.39) 

ro 

Here we assume that A''i = Ni{x,t), N2 = N2{x,t), but that iVg = Nsit), i.e. the boundary 
inhomogeneity only depends on time. In order to use the linear theory we have developed, 
we must rewrite this as system with homogeneous boundary conditions. To accomplish this 
we will utilize the following lemma. 

Lemma 3.4. Let 

Hx,t) = -^rl{x). (3.40) 

Then for each t, ip(t) satisfies C2,^{t) = for x E (0, M) and B^{t) = N^it) at x = M. 
Also, Ciiljit) = -NBit)poix)/5. 

Proof. The results follow from simple computations. □ 

With this ip in hand, we can reformulate fl3.38p - fl3.39p so that the resulting problem has 
homogeneous boundary conditions. Let w = ip + w. Then Lemma [3.41 implies that 

9. (1) ^ c (:^ J^A. (z^i) ^ ( J; t ^^s'L^ (3.41) 



w / \w 



N2J \ J \N2 + dtNsrl/{36) 



along with the boundary conditions 

w 



(0, t) = a{M, t) = and B{w) =0 at x = M. (3.42) 
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Employing the variation of parameters, we can then solve fl3.4ip - fl3.42p via 



ait) 
w{t) 



..tc 



a(0) 
w{0) 



Nlis] 



ds + - 





rlMfYh,)ds- (3.43) 
dtNB{s)r^/3j 



We can then go back to w = ip + w: 



a{t) 



AC 



a(0)\_l/ 



+ 



N2{s] 



ds + - 





(t-s)c ( Nb{s)pq 
dtNB{s)rl/?> 



ds. (3.44) 



Now let us define a norm for the pair a, w given by 



M 



a 
Po 



dx -\ — 
2 



M 



W\ 



-dx 



^\ la ^^^^ {^'^^ \9xw\^ + ypo rld^ 



We also define 







M 







Po 



dx 



M 



\dtw\' 



■dx. 



dx. (3.45) 



(3.46) 



We can then recast the result of Theorem 13.31 as 

2 



o-(O) 
w{0) 



< Ce'^'(B{a{0),w{0)). 



(3.47) 







Using these quantities and estimate f l3.47p . we can record estimates for solutions to f l3.38p - 

(Km . 

Theorem 3.5. Suppose that a andw solve the inhomogeneous linear problem f l3.38p - fl3.39p . 
Let ip be given by Lemma \3.4\ and w = w — ip . Let ||-||q and £(■, ■) be given by fl3.45p and 
f l3.46p . respectively. Then 



ait) 
wit) 



— e 



tc 



< / Ce^^'-'W^iNiis),N2is))d.s 



aiO) 
wiO) 

+ J \NBit)\ + T fe'^'~'^ {\Nsis)\ + \dtNsis)\ + \d^Nsis)\) ds. (3.48) 



6 



Proof. From the above analysis, we know that a and w are given by (13.441) . where e*^ is the 
homogeneous solution operator given by (13.371) . Hence fl3.47p implies that 



ait) 
wit) 



AC 



aiO) 
wiO) 



1 



< 

-'o 


Nsit)rl/3 



Ce^(*-^) ,/(BiN^is),N2is))d.s 



^\fo \j<^iNBis)p,,dtNBis)rl/?>)ds. (3.49) 
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Then, since N^lt) is only a function of time, not of x, we can easily estimate 



<C\NBit)\ 



(3.50) 



and 



(E{NB{s)po,dtNB{s)rl/3) < C {\Ns{s)\ + \dtNB{s)\ + \d^Ns{s)\) , (3.51) 



where C > in fl3.50p - fl3.5ip is a constant depending on various (finite) integrals of po and 
tq. The estimate f ]3.48p then follows by combining f l3.49p -f !3.5ip . 

□ 



4 Nonlinear energy estimates 
4.1 Definitions 

We are interested in small perturbations a, v around the stationary solution p = po, r = rQ, 
and V = 0. In particular, we assume that 



9 11 



(4.1) 



This assumption will be justified later when we close the nonlinear energy estimates. For 
such small solutions, the Navier-Stokes-Poisson system f ll.l7p and fll.lSp can be written as 
follows: 



dta + Anp'd^ir'v) = 
dtv + Anr^d^P - Anrld^^Po + ^ - ^ = 167T^r^d.,iiAe/3 + 5)pd.,ir^v)). ^^'^^ 



The dynamics of r are determined by 



dy 



1/3 



and dtr{x,t) = v{x,t). 



(4.3) 



It turns out that it is convenient to analyze a/po rather than a itself, so we rewrite the 
continuity equation as 

^d,(-]+ATrpd^ir%) = 0. (4.4) 



(4.5) 



P \Po, 

We will also rewrite the momentum equation. To do so, we first note that 

Anr^d^P - Anrld^Po + ^ - ^ = Anr'd^iP - P,) + x ( ^ - ^ 



and then note that for small perturbations satisfying (14. ip . P — Pq = K{p^ — Pq) can be 
written as 



P-P^ = pI{ Kj- + a 
Po 



a 



Po 



(4.6) 
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where a* is the smooth bounded remainder from the Taylor's Theorem. We then rewrite the 
momentum equation as 



a 



dtv + Airr { K^pl h a^pl 

Po 



1 



V, 



(4.7) 



where V := 167r^r^93;((4e/3 + S)pdx{r'^v)). We give an equivalent expression for V so that it 
appreciates the boundary condition (11.201) in energy estimates: 



where 



3 \r 



(4.8) 
(4.9) 



satisfies W(M) = because of the boundary condition (ll.20p . We use v to denote the 
minimal viscosity coefficient: 

v := min 15,— 



(4.10) 

We now define instant energy functionals for a and v. In what follows, all of the integrals 
are understood to be over the interval [0, M]. 



8^ 



8^ 
8' 
8^ 



\v?' dx -\ — 



'1 + ^)2 

po' 



a 
Po 



da; + - / u 



To 



dx ■= 8°'" + 8°'" + 8^^'' 



S j 167rV|5,(r\)|'t^a; + y j lOvrV^ 





2 




dx 







6 + —] 167rV 



.2 4 



1 



2 1 

\dtv\ dx + 



1 + - 

po 

7-1 



Kipl 



2J (1 + ^)2 



Po 
a 
Po 



dx := 8^'" + 8^''' 



dt 



dx := ^2'- + ^2,<x 



5 j lQ'K'^p\d^{r^dtv)\^ dx + ^ j 16 



n'^pr^ 





2 

dx 







(5 + ^1 4npo 



a 



dx{r d^[ — 
\Po 



dx 
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The corresponding dissipations are given by 



dx 



po 



:= j \dtv\'dx + j Wn^K-fr^pl 
p3 := j \d^v\^dx 



dx := V^'"" + V^'" 

2 

dx 



2 ^^7 



a 



dx{r d^[ — 



dx. 



In addition, we introduce various bootstrapped and auxiliary energies and dissipations (de- 
noted with subscript b and a, respectively) that can be controlled with the above instant 
energies and dissipations: 



0,r 



r 



dx 



8, 



1,(T 



5 + —1 167r 



3 / 



cr 



Po 



dx 



(4.11) 



l,cr 



2^7-1 



03,0- 



levr^ 



dxdt 



£. 



3,1; 



— \dx{pdx{r'^v))\^dx] 8, 



a 
Po 

3,1- 

112 



Po 



dx 



dx 



pr 



(4.12) 



dx 



8^ - 



5 + ^1 47rpo 



{ r^dtdx { — 



dx 



8t^-= j ISttVo \ dx{r'^dx{pdx{r'^v)))\^ dx. 
Finally, we introduce some bootstrap energies that depend on a parameter (5 G 



(4.13) 



8^ 



0,0- 



/pr 


o 


/ p 


po 



p 



Po 



dx. 



(4.14) 



For the proof of our instability in Section O we will need to invoke higher order energy 
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functionals and dissipations, which we define now. For i = 2 and 3, let 

2 



PO' 

i2 , 4e 



dx 



6 I lQ'K^p\d^{r^div)\ dx+ — / IGttV^ 



5 / 167rV|<9,(r29^t;)|'rfx + y / IGttV^ 



I 9^+^17 1 (ix. 



r 



Next we define bootstrapped energies and auxihary energies for i = 2 and 3 



(4.15) 



1 

Po 



Po 



l+2j,t) ._ 



6 + j\ levrV 



c2+2i,a ._ 

The then define the total energy by 



Po 



(ix. 



(4.16) 



i=0 



(4.17) 



Throughout the rest of the section, we assume that 



a 
Po 



dt 



a 



Po 



a 



Po 



r 



pr'd, - 
Po 



a 



+ 



d 



a 



Po 



r 



dt'v 



dtv 



< Oi (4.1^ 



for sufficiently small constant di, where the norm IHI^^oo is over the spatial region [0,M]. 
The validity of this assumption within the total energy S will be justified in Lemma 14.91 

Since r is determined through an integral of a as in (14. 3 p , for small perturbations satisfying 
(14. ip we may use Taylor's Theorem to write — as 



'"o 



1 + 



—dy+ — 



'^o ^0 P* VPo 



cr 



dy + 



C2 



^dy 
Po 



(4.19) 



where p*/po ~ 1 is a bounded smooth function of ^. Hence the 1 — tq/t estimate (up to a 
constant) in (I4.18P can actually be guaranteed by the smallness of the other terms in (I4.18p . 

The relation (I4.19p will be useful in various places. We now record a couple other useful 
identities. 
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Dynamics of y : From (14. 3 P we have that 

'rQ\ VqV fv\ /ro 



a 



r 
r 



1 



1 



'ro\3 



+ 



PO 



Some useful inequalities and identities: For any v (not just solutions), 

2 

pr 



V 4:71 

r 3 



^pda:{r'^v) - pr^d^ (^^^ 



p \ dx{r'^v)\'^ + pr^ 



dx 



pdxirv"^) = pd^ 



2-pdx{r v)-- 

r Air r-^ 



V ( / V 

- ipdxir^v) + pr^dx - 
r I \r 



(4.20) 



(4.21) 



4.2 Estimates 

Throughout the rest of the section, we use C to denote a generic constant that may differ 
from hne to hne, and t] to denote a sufficiently small fixed constant which will be determined 
later. The constants C are allowed to depend on rj, which presents no trouble in our ultimate 
analysis since first we will fix an rj, which then fixes the constants. 

In the following series of lemmas, we provide the energy inequalities for S. We start with 
£° and V^. 



Lemma 4.1. 



;^^° + i^°<C(i + ^i)£:° + V 



Proof. Multiply (14. 7p by v and integrate to get 



(4.22) 



--r I Ivl"^ dx ■ 

2dt 



ATTdxirh) { K-ipl— + ( — ) \dx+ I v 



Po 



a 



Po 



For [i] we use (14. 4p to see that 
Po 



-A 



a 



Po 



dx 



1 d 
2Jt 



K-fp 



,7-1 




^ PO ' 



a 
Po 



dx + 



(7^7 + (1 + ^)a.)pl 
(1 + -)' 



,7-1 



x(rQ — r^) 



2ri 



dx 



j vVdx. (4.23) 



(Hi) 



a 


2 






9,1 




Po 







(4.24) 
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However, 



{K^ + (1 + ^)a.)pg 

(1 + ^)3 



,7-1 



cr 
Po 



dt [ — ] dx 

Po 



< C(l + 



(4.25) 



For (m), the Cauchy-Schwarz inequality yields 

|2 



J pr^ V ' 

9 



X 



-^{r +ro)(r + ro) 







r 



(4.26) 



where we have used (14.211) at the second inequality. From (14. 8 p and the boundary condition 
W(M, t) = 0, we get 



(m) 



-5 J 16n'^p\d^{r%)\^dx-^ J Wn^pr^ 0^ 



dx = -V^ 



(4.27) 



Next, from (020]) 

J^d_ /" -| _ ^0 
2dt 



vro 



dx = —V I 1 ] dx 

\v 



<v I ~^d,x + u I p 



r 



r 



dx < + C£^'\ (4.2^ 



The desired estimate then follows by combining these estimates. 



□ 



With Lemma 14.11 in hand, we can bootstrap to control ^ with an improved weight. 



Multiply (14. 4 p by p^-^ and integrate to get 

Po ^ 



f P^—dt (—\ dx = - [ pho— ■ 4:7Tphjr^v)dx. (4.29) 
J P Po \PoJ J Po 



Thus 



Id fP^ ^\^^C /■ 16^2^ 1 (^2^) 1 2^^^^ /■ 2/3 

J P Po V J J 



2dtJ p 



dx 



a 

Po 



1 fpr 



2 J p' 



dto- 



o 

Po 



dx, (4.30) 



which means that 



(4.31) 



Next we consider and V^. 
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Lemma 4.2. We have that 



2,v 



where 



q:=qi + qi 

is bounded due to fl4.18p . 

Proof. We divide the proof into steps. 
Step 1 - S^'" and V^'" 
Multiply (14. 7p by dtv and integrate to get 



1 + - 

PO 



dx 



{iv) 



+ f dtv'^^^^dx- [ dtvVdx 
J r^r^ J 



(v) 



(vi) 



For (iv) we first expand 

(tv) = J Anr^dtv S^K^dM)y^ + K^pld, (^j^ | dx 

+ j Anr'^dtv I (a^Po) (^y^ + 2a^plyd^ (^y^ I dx 



and then estimate 



[iv)i + iiv)2 + (if Is + (iv) 



1 



M, + {iv)3\ < - I \dtv\^dx + c{i + el) 



/ 


a 




Po 



dx < -V^''' + 
4 



and 



\{iv)2 + {iv)i\ < / M'K^r^'pl 



Po 



-a' (- 



dx 

2^ 



pI \dtvfdx 
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For (v), we get 



{v)\ < / \dtv\''dx + 



X 



(r +ro)(r + ro) 



|r-ro|'rfx < ip^''' + C£°'^ (4.3^ 



The term (vi) forms the energy £ ''" and nonhnear commutators: 
Id 



1 

+ 2 



+ 167r29^(f • 2rv)pd^{r^v)dx+^ J 167i^pr^d^ (v ■ (-^)) (^) t^a;. 
Using f l4.18p and the fact that 

d,{v ■ 2rv) = 2 (^d,{r\) + rhd,{^)) 



(4.39) 



dr. 



(4.40) 



the absolute values of the second and third lines may be bounded by C9i£^'^. 
We may now combine the above to deduce that 



dt - 2 



(4.41) 



Step 2 - S^'" and V^''' 
For the estimate of dx{-^^), we first rewrite (14. 7p by replacing pdx{r'^v) in V by f^t(^, 
through the continuity equation (14.41) : 



4£ 



P \Po 
a 



S+-\ Anr' ^dA {-]+d,{^\dt{-\} + dtV + 



Po 



a 



Po 



x{rQ — r"^) 



Po 



+ K^d, {p2) - + d^ia^pl) (-] + 2a,p^-d, ( -] \ = 0. 
Po \Po/ Po \Po/ j 

(4.42) 



Note that d^if) = -{1 + fj-^d^i-^). Multiplying by Airr'^d^i-^) and integrating, we 

are led to the estimate 



1 d 



1 + - 

PO 



Po 



\Po 



dx + j 167i^K-fr^p2 dx ' 



dx 



(4.43) 



2r^{5 + f 



"jrAdM'dx. 



Note that the last term in ( 14.43^ may be bounded by q2^'^'^- We then obtain (14.32^ by 
combining (I4.4ip and (I4.43p . □ 
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The estimate (I4.32p is not of a closed form by itself. Its use will come by coupling to the 
result of the following lemma. 



Lemma 4.3. 

d 



(Jjlj ^ 



(4.44) 



Proof. We take dt of (|12D to see that 



d^v + 4nr'dAKjp2dt 



+8TTrv < K'yp^dx 



a 



Po 



a 



Po 



a 



2a^pl—dt ( — ) + dta^pl 



Po \Po 

2 

a 



Kl9x (Po) — + dx{a^pl 

Po \Po 



a 



2a^pl—d^ 



a 



Po \Po 
xf (rg — r^) ^xrv 



(4.45) 



The energy estimate fl4.44p may be derived from 04.451) as in Lemma 14.11 we multiply fl4.45p 
by dtv, integrate over x G [0,M], and integrate various terms by parts in order to identify 
dS'^/dt, T>^, and some error (lower order or commutator) terms, the latter of which may be 
estimated by the right side of fl4.44p . Since the argument is essentially the same as that of 
Lemma 14.11 we present only a sketch. 

First observe that the product of dtv with the first two terms in the first line in fl4.45p 
forms the energy term dtS^ and some error terms: 



dtV 



Po 



dx 



PO' 



\P0 



dx} + Z, (4.46) 



where 2^ is a term whose absolute value may be estimated by the right side of (I4.44p . Here 
we have used the continuity equation fl4.4p and an integration by parts on the second term. 
Next, we compute 



dtV = leTrV^a^^tW + 16n\2rv)d^W + —UTrr^d^ — 

3 \ r 



(4.47) 



and note that the boundary condition W(M, t) = implies that dtW{M, t) = as well. This 
allows us to integrate by parts without introducing boundary terms: 



IQn'^r'^dxdtWdtvdx 



lQ7T'{r'dtv)dtWdx. 



Using this, we find that 

j dtVdtvdx = -5 J 167rV|<9^(r^<9tt;)|2rfx+ y J 16n^pr' 
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a. I*; 

r 



(4.48) 



dx + Z, (4.49) 



where again Z is an error term with the property that \Z\ is bounded by the right side of 

Finally, all of the remaining terms that arise when we multiply fl4.45p by dtv can also be 
estimated by the right side of (14 .44 p . For example, the second term in the third line can be 
estimated by noting that xr jr^ is bounded, which means that 



- y A^dtvdx <7] j IdM'^dx + C j |t;pda; < 7]£^ + C£^. 
Combining all of this, we find that 



(4.50) 



PO' 



Po 



dx 



r 



(4.51) 



which yields (I4.44p . 

We now derive bootstrapped estimates for dt{-^)- We take dt of (14.40 to get 



pI 



a 



P \PoJ P^ \ \Po 

Next, we multiply (I4.52p by Po^t(^) and integrate to see that 

Id fp^^' 



□ 



(4.52) 



cr 



Po 



dx < 



C 



IQn'^ p\dx{r'^dtv)^ dx + 7] / pp^^^ 



a 



Po 



dx 

3 



2dtJ p 

- + Pr'a. (^) } 43. (^) + I / (a, (^) ) 

-An j dtad,{r^v)pldt (^j^ dx. (4.53) 



Then we estimate 



j Svr^ {pd^{r\) + pr^d^ (^) } p^9, (-^^ dx 

-ATI j dtad,{r\)pldt (^j^ dx < C9iV' + C9i£^f^, (4.54) 



to obtain 



dt - A 
Next we estimate £^ and V^. 



(4.55) 
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Lemma 4.4. There exists an energy J-"^ so that 



d 



dV 



3 e, 



8 4 



'1 I 



Moreover, we have the estimate \J^^\ < C6i{£^ + 
Proof. First recall fl4.45p and rewrite it as 



dx{a^pl)( — ] +2{a^pl)—da^[ — ] 
\PoJ Po \PoJ 



+ STirv { K^pld,{-) + K^dM)- + d,{a,pl){-f + 2a,pl-d,{-] 



Po 



Po 



Po 



po Po 



(4.56) 



(4.57) 



where dtV is given in fl4.47p . To derive fl4.56p we will multiply by dfv and integrate over x. 
We divide the estimates into the following steps. 
Step 1 - (ai), (02), (6), and (c) 

We begin with an estimate of the product of d^v with the terms (ai), (02), (b), and (c). 
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First, we use fl4.42p to replace dtd^i-^) by lower order terms: 



{Pl)dt 



+ 47rr^ < dt 



dx{.a^pl)( — ] +2{a^pl)—d^(—] 
\PoJ Po \poJ 



a 



P 



PoJ Po 



K7 + 2a,-)pg f 

PoJ po{S+f) 



dt 



dtv + ^ ^ , + 47rr2 



Klppx { — 
\Po 



+K'yd^{p]^)— + d:^{a^pl) { —] + 2a,p2^— "^x ( — ) > 
Po \PoJ Po \PoJ\) 



dx{a^Po) 



a 



+ 2pldt a*— ]d^\ — 



Po 



a 



\Po 



iA^) + (A,) + (A3). 



Then J d^v ■ [(oi) + {a2)]dx can be estimated as follows: 



(4.58) 



d^v ■ {Ai)dx < 



dfvfdx + COi 



Po ^4 

P' 



dx < + COiS 



1,0- 



J d^v ■ iA2)dx \d^v\^dx + C[£^''' + + (1 + ^i)(^''" + ^°'")] 

J dfv ■ {A:^)dx \dtv\^dx + CdiE^'" + C{1 + ^i)^o'". 

Next for (6) and (c), we may estimate 



(4.59) 



d'tv ■{h)dx<- / \d'tv\'dx + C^i(£i'" + 



y a^^; .{c)dx<^ J \d^v\^dx + C(£°'" + S"^''). 

Combining the above, we arrive at an estimate for J dfv ■ [(ai) + (02) + (&) + {c)]dx. 
Step 2 - The viscosity term 

Now we consider the viscosity term, dtV. We claim that there exist J^^, G so that 



(4.60) 



d 



d 



d^v ■ dtVdx = -—£'^ - — ^ Q 



dt 



dt 



where 



< Cei{S^ + S^) and 
\G\ < ^V'' + Cei{£'' + £l''' + £^). 



(4.61) 



(4.62) 



38 



Recall that dtV may be computed as in f l4.47p . and that dt}V{M,t) = 0. Then a simple 
but lengthy computation, using integration by parts, reveals that 



dx 



+ Go + Y, 



where 



(4.63) 



Go 



2p r 



4:6 _ 



dtV 



dx 



+ jpr^d. {^) [s32Air^0,v) - f IGvr^^r^^ (^) } dx, 



(4.64) 



and Y = Yi + Y2 with 



= -iQn^ j [5dt(Td.^{r'^v) + 5pd^{2rv'^)\ d^{r^d^v)dx 



lQn^J\%,ipr')d.{^)-^pr'dJ^ 



Y2 = 327r2 / rvd'}vd^ 



6pd^{r^v) + ^pr^S^ (- 
3 Vr 



d^{r^d^v)dx, 
dx. 



(4.65) 



Let us define the quantity Q so that Yi = —IGn'^ f Qdx{r'^d^v)dx, i.e. Q is the sum of the 
bracketed terms in the Yi integrand. Then we may compute 



1^1 = ^ / -16TT^d^{r^dtv)Qdx + I 167r^d^i2rvdtv)Q + d^,ir^dtv)dtQ) dx 



(4.66) 



dt 



Similarly, we have that 

Y2 = j^ f -167r'd^{r'd,v)y 



4:S / V 

5pd^{r'^v) + —pr^d^ - 
3 Vr 



dx 



2v /2?A 2v 

dJ2rvdtv)—W + d^{r^dtv)dt — W + d^{r^dtv)—dtW 
r \ r J r 



dx 



(?) 



(4.67) 



- 167r^ j r^d'fvWd^ [ — ] dx 
d 

di' 



-^J'l + G^. 



Combining the above, we find that (14.611) holds with J-"^ = J-'f + J-f and G = Gq + Gi + 
G2- To complete the proof of the claim, we note that the estimates (14.621) follow from the 
definition of J-"^ and G, using (14.52^ to replace d^a by other terms. 

Step 3 - Conclusion 
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The only term that remains is 



j dlvdlvdx = 



(4.68) 



With this, all of the terms in fl4.57p are accounted for. We may then combine the analysis 
of Steps 1 and 2 to deduce the estimate fl4.56p . 

□ 



We now bootstrap more estimates. First, we multiply ( I4.20p by ^(1 — ^) and integrate 



to get 



(4.69) 



where we have used (14.2 ip to control J \vf /{r^p)dx < CS^'"". Next, by multiplying the 
equation (14.42P by ^f^s(^) and integrating, we get 



dt 



(4.70) 



Note that here we have again used (I4.2ip to control / \dtvf jir^p^dx, which is possible since 
(I4.2ip is valid for any choice of not just solutions. From the equation (I4.42p we also see 
that 



3,cr 



3 / 



P 



drr di ( — 

Po 



dx < C{£-' + El''' + + 



(4.71) 



Next, by applying dx to equation (14. 4p . we find that 



P \Po 
We then use this to get 



" po 



+ A'Kdx {pdx{r\)) = 0. 



(4.72) 



= / - IdApdAr'v)) f dx < C{S!''' + 8^- 
Since / pr^\dx{^)\'^dx < C£^, the equation (14.70 implies that 

= J pr^l^, (pr^S, (^))|'rfx<C(^2 + £l+0• 



(4.73) 



(4.74) 



We now illustrate how the higher order energy estimates of spatial derivatives of dx{ — 



and dx{pdx{r'^v)) work. The following lemma concerns the estimate of dx{r^dx{-j^))- 



PO' 



Lemma 4.5. 

d 



+ < (r/ + C^O^ ^ + C(f ^ + + £l'' + + Cei(^,''" + £l 



dt 



'1,0" I 473, (T\ 



(4.75) 
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Proof. First, we multiply fl4.42p by and apply to get 



"3 



P 



(i) 



(n) 



P 



a 



Po 



(^0 - r^) 



iiii) 



(iv) 



xr 



Po 



(v) 



+ 



(vii) 



[ Po VPo/ VPo/ 

+ 47r ^2a^pljd,, [r'd,, (^j^^ + 2r^9. («*/^o^) (;^) } = 0' (^-^S) 

The energy inequality fl4.75p can be derived as in Step 2 of Lemma Wf2\ by multiplying f l4.76p 
by pdx{r^dx{-^)) and integrating over x. We provide the details on how {i) — {vii) can be 
treated; other terms can be estimated similarly. 



dx { r'^dx - 
VPo 



dx 



S+j^AttJ podx {dt{r')dx ) dx {r'dx [j^ ) dx . (4.77) 



(*) 



Since dt{r^) = A^r\ 



4-po 



/ \ \ 2 

4r. / O- 



dx{r dx[ — 
VPo 



dx + J Adx (^) r'dx (^^^ podx (^r% (^^^^ dx, (4.78) 



and since - and pr'^9a;(-) are bounded by 9i, 



For (ii) we write 



[ii) ■ pdx [ r^dx \ y \ \ dx 



a 



*)\<ce^{s^ + s. 



1,CT\ 



(4.79) 



1 + 



po 



—=dtdx 
y/P VPo 



y/pdx r a 



po 



(4.80) 
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and therefore, 



It is easy to see that 



cr 



[ii) ■ pd^{r d^{ — ))dx 



Po 



(4.81) 



a 



[{in) + (iv) + (v) - (vi)] ■ pd^ [j^ ) ) 



Finally, (vii) forms the dissipation V^. 

We also get an estimate for dx{r^dtdx{-^)) from fl4.76p : 



(4.82) 



□ 



5 + y ) 47rpo 



dx ( r^dtdx ( — 



dx 



(4.83) 



To derive an estimate of the third spatial derivatives of f , we first multiply fl4.72p by 
and then apply dx'- 



'-dx [ r'dxdt (-]] +2T'dx ['^\dxdA-\+dx[T'dx['-\\dt[- 



Po, 

p 



Thus, we obtain 



VPo 



f \ . f Ar. f POW f O- 

Vpo/ V VP// Vpo 

+A7idx {r^dx {pdx{r\))) = 0. 



(4.84) 



S. 



0,2 



167r% \dx{r''dx{pdx{r'v)))Y dx < El + B\{EY + E^). 



(4.85) 



We now present the higher order energy estimates. We start with E'^ and E\. 
Lemma 4.6. 



+ C^i(^f'" + + £2 + £^) + 



(4.86) 
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Proof. We apply dt to fl4.45p to see that 



cr 



Po 



a 



Po 



Pi 



Po 



cr 



Po 



+Adta,-dt — +dta, — 



Po 



+ 167Tvrdt 



( o" 



'9x(a*Po)(— ) +2{a^pl)—d^[—\ 
\PoJ Po VPo/ 



VPo 



cr 



/ cr 



Po 



VPo 



cr 



Po VPo 



2x{rQ — r'*)9tf 4:xrdtv 12xf ^ ^ 6a;(rQ — r^)v^ g^-j 



4^4 



We derive the energy estimate of f l4.86p from f l4.45p by proceeding as in the proofs of Lem- 
mas 14.11 and 14.31 That is, we multiply the resulting equation by d^v and integrate over x, 
integrating by parts in some terms to recover dS^/dt, V^, and some error terms that can 
be estimated by the right side of f l4.86p . Since the method of proof is already recorded in 
Lemmas 14.11 and 14.31 we omit further details. 



□ 



An estimate of dx{pdx{dt[r^v])) can be obtained through the equation (14.451) : 
We now bootstrap to control d^{-^). We apply dt to f l4.52p to get 



(4i 



VPo 



An—pd^{r^d'fv) - 2Ati— pd^{rvdtv) -Sn— pd^{v^) 

(a) (b) 



d+cr ^r, ( o 
+ 6—92 — 

P VPo 



6 



PoP^ 



(4.89) 



Note that (a) = Ipd^ir'^dtv) + pT^d,{^^-f and (6) = 'ir^pirfQn) + X(H)3 and thus by 



multiplying (I4.89P by ^"^^l^) and integrating, we obtain 
dt 



(4.90) 
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Next, we take dt of (14.421) to see that 



5 + %]Airr^\P^dld^ 



3 



P 



Po 



dtd^ 



dt 



^o_^2 fPo^^ f 



3 / r 



P \Po 



a 
Po 
Po 
P 



V X xr 



+ dtv-2-{- + — + Anr' K^pldtd. - + K^dM)dt - 



Po 



Po 



\PoJ Po \PoJ 

r [ \PoJ Po \PoJ Po \Po 



= 0. 
(4.91) 



Therefore, by squaring (14.910 and integrating, we find that 



5 + l£) i6vrV 



\Po 



dx < C{S^ + S') + Ce\{ej[ + El'' + f 1 + ^o'")- (4-92) 



Also, by first dividing (I4.9ip by r and then squaring, we obtain 



— \ IGttV 



a 



\Po 



dx < Ci£^ + S'l'' + E^jl + 



(4.93) 



Now we record an estimate of £^ . 
Lemma 4.7. There exists an so that 

j^[s^+F] + < + ce^)£' + (i + ce,)v^ 

Moreover, < Cei{£^ + £^ + £^). 



(4.94) 



Proof. The energy inequahty (14.941) can be obtained by multiplying (I4.87P by dfv and inte- 
grating over X as done in Lemma 14. 4[ We omit further details. □ 

As seen in the previous estimates in Lemmas 14. 31 14.61 14.41 and 14. 71 the time differentiation 
of the equation keeps the main structure of the highest order terms as well as the boundary 
condition. Using the time differentiated equations (I4.87P and (I4.89p . we can follow the line 
of analysis presented in these four lemmas to derive energy inequalities for 8"^, Si, £^I'f , £^ 
and Sa'''- We record these in the following lemma but omit a proof. 
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Lemma 4.8. Let £ he given by (14.171) . We have the following estimates. 

+ COi^Sl^'' + 8^ + + 8"^ + 8^ + 8^) + C8^ 
8l < C{8' + 8l'^ + 8^) + Ce,{8l'^ + 8l'^ + 8^ + 8^ + 8l'^ + 8l) 



dt 

d 



1 



— [8^ + J^«] + I)« < (r^ + C0i)£« + (- + + C{8 - - C") + 

(JjL ^ 

where \J^^\ < 09^8 + C\8\^ 



(4.95) 

(4.96) 
(4.97) 

(4.98) 
(4.99) 



The next lemma ensures that the assumption (I4.18P is vahd within our energy 8. 
Lemma 4.9. There exists a constant n > so that if 8 < k, then 



a 



+ 



+ 



1 - 



+ 



pr'd. (- 



+ 



a?(- 

\P0 
+ 



dtv 



+ 



dtV 



< CVS (4.100) 



for some constant C > 0. Here 8 is given by (I4.17p . 



Proof. The proof proceeds through several steps. 
Step 1 - d^{a/po) estimates 

We begin by estimating a/po in W^'^{{0, M)). First, we use Holder's inequality to esti- 
mate 





a 


dx < Vm (^j 


o 


dx^ 


/ 








Po 




Po 





1/2 



(4.101) 



On the other hand, we may estimate 



/ 7 r (;^) I ^7 (/ ^^') - ^ ^^''^^''^ - ^'-'''^ 

Here we have used the fact that r^(x) ~ x^/^ for x ~ 0, which follows from the definition of 
r(x) and L'Hospital's rule, to see that / {p/r'^)dx < oo. Combining these estimates with the 
usual 1 - D Sobolev embedding W^'\{0,M)) ^ C^{{0,M)), we find that a/po E C° and 



a 
Po 



(4.103) 



L°° 



Now to control dtc/po we argue similarly to estimate 

dec (^y-^ dx < C^/8^ + < CV8. (4.104) 
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/ 


dta 


+ 






Po 







Then dtO j E C° and 



dta 



Po 



< cVs. 



(4.105) 



A similar argument, employing and for i = 1,2, then implies that 









+ 


Po 




Po 



d^a/po,dfa/po E C° and 



We thus deduce from (141031) and flirTOSD - flCTe]) that 



(4.106) 



a 

Po 



+ 



cr 



\Po 



+ 



L°° 



+ 



a 
\Po 



< cVs 



(4.107) 



Step 2 - 1 — ro/r estimate 

Let us now suppose that £ < k with k small enough so that C\/S < |, where C > is the 
constant appearing on the right side of (14.1071) . In particular, this implies that ||o"/po|Iloo < 
1/2 < 1. With this estimate in hand, we can derive an estimate for rg/r. Indeed, the Taylor 
expansion (14.19!) easily implies the estimate 



r 



< C 



a 
Po 



l+k 



< cVs 



(4.108) 



for some A; > 0. This is the 1 — tq/t estimate in (I4.100p . 
Step 3 - d^v/r estimates 

We now turn to estimates for d^vlr, = 0, 1, 2. From Step 1, we know that — and df — 

til 11 f 1 pQ 5pg 

are continuous and bounded, while from Step 2 we know that — < 1/2 so that 1 + — 

is also continuous and bounded. From the boundary conditions at x = 0, we also have that 
r^T;(0,t) = 0. Hence we may spatially integrate the continuity equation (14. 4p to see that 



-1 



(r v){x, t) = — 



An 



Notice now that, due to the asymptotics (ll.24p . we have that 

dy 



dy. 



(4.109) 



< oo. 



'o Poiy) 

This and the estimates (14.107!) then imply that v/r E C° and 

< cVs. 



(4.110) 



< C 



dt 



a 



Po 



(4.111) 



On the other hand, due to L'Hospital, we have that 



1 



^' dy 4:7rp{x,t) Air 



r^ix,t) Jo poiv) Spoix) 



^ 1 + 



a{x, t) 
Po{y) 



< cxo for X ~ 0, 



(4.112) 
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which means that 



1 



xG{o,Af) r^(x,t) Jo Po[y) 



^ dy 



< oo. 



We may then deduce that v/r G and 



V 


< C 






9,- 


r 




Po 



1 



dy 



sup ^r—. I t X 
LOO xe(o,M) r^(x,t) Jo Po(Z/) 



Now we apply dt to f l4.4p and argue as above to see that 

1 d^y^t) 



(4.113) 



(4.114) 



(r dtv){x,t) 



+ 



dt(riy,t) 



Poiv) 



-dy 



dy 







'^2(r2t;)(y,t)4^rf2/. (4.115) 



Using this, we may argue as above (using the estimates (I4.11ip and (14.114^ ) to deduce that 
r'^dtv,dtv/r G C° and 



dtv 



< cVs. 



(4.116) 



An iterative argument, using apphed to (14.41) in conjunction with the estimates (I4.116p . 
then allows us to see that r'^d'^v, d^v/r G C° with 



(4.117) 



Then fl4.11ip . f l4.114p . and (I4.116p - (l4.117p may be combined to derive the d^v/r estimates 
recorded in (I4.100p . 

Step 4 - pr^dx{<7/ po) estimate 

Since ||a"/po|lLoo < 1/2, to prove the pr^dx{cr/ po) estimate listed in f l4.100p . it suffices to 
estimate this term with p replaced by po- We claim that 



Por'dxi-) 




Po 





(4.118) 



To prove fl4.118p . we will use the 1 — D Sobolev embedding W^'^ ^ First note that 



Por 



dx 



dx < 



pplr'^dx 



1/2 



r 
P 



dx 



1/2 



dx 



On the other hand, we may compute 



dx Por dx — 



a 



Po 



.Po 
and dxPo 



a 



Po 



'-dx [r'dx[-]]+ dxPor'dx ^- 



a 



'1,0- 



(4.119) 



47rp VPo 



X 



(4.120) 
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Then since -^St=^ < ^ as long as 7 < 2, we may estimate 



dx Por^dx 



a 



Po 



dx 



< C / porl 



Pi 



a 



Po 



2 2 
P' 



a 



Po 



Then from this and Holder's inequality we get 

dx ^ 



dx por 



a 



Po 



dx < 



Porl 



Porl 



2 2 6 

X r 

+ 



27-2 8 

Po "^0 



dx\ 



( 



\P0 



dx 



<CE^^El\ (4.121^ 



dx por^dx 



a 



Po 




<CVS'^ + Jsl'\ (4.122) 



Together, the estimates fl4.120p and f l4.122p constitute a W^'^ estimate for Por^dx{cr/ pq), so 
we then obtain fl4.118p via the Sobolev embedding. □ 



5 Nonlinear instability 
5.1 The bootstrap argument 

Based on the nonlinear estimates in the previous section, we now establish a bootstrap 
argument that allows us to control the growth of S in terms of the linear growth rate A, 
constructed in Theorem 12.11 The idea is to assume small data and that the lowest order 
energy, S^, grows no faster than the linear growth rate; then the inequalities in the last 
section allow for a bootstrap argument that shows that all of £ grows no faster than the 
linear growth rate. 

Proposition 5.1. Let a and v be a solution to the Navier-Stokes-Poisson system (14. 2p . 
Assume that 

v^(0) < Cql and \^{t) < Coie^^ for < t < T, (5.1) 

where £^ and £ are as defined in (14. ip and (I4.17p . Then there exist and 6'^ > such that 
ifO<t< min{T, T{l, 9^)}, then 



V£{t) < de^' < CJ., 



(5.2) 



where we have written T{l, 6^, 



Proof. To prove the result we will employ a bootstrap argument using all of the nonlinear 
energy estimates derived in the previous section. We now choose 61 and rj sufficiently small 
in all of these estimates so that C9i + r] < ^ and C9i < 1/8 in all of the energy inequalities. 
Throughout this proof we will write C for a generic constant; we write this in place of C to 
distinguish the constants from those appearing in the nonlinear energy estimates. 
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To begin the bootstrapping, we show that the estimate (15 .ip allows us to control an inte- 
gral of the dissipation. Indeed, we use fl4.22p and flS.ip along with Gronwall's inequality 
to see that for < t < T, 

^go ^ ipo < (j^o < cC^Me^'y + -£'^- f e''/^^'-'^V%s)ds < Ci^e^^K (5.3) 
dt 2 - - ' 2 2 Jo V ; - V ; 

Then we employ (I4.3ip with /3 = — 1 in conjunction with (15. 3p to see that 

(JjL 

Jo 

^!'r(t) < Ci^e^^' S^'^'it) < C^^f (t) < Ct^e^^'. (5.4) 

Next, let g > be the constant from estimate (14.321) and choose k = (4g)/A. Then (14.321) 
and (I4.44p . together with (I4.22p and the above estimates, imply that < t < T 

^[ks' + ks' + s'] + < + cei){s' + ks^) + qs^ + c{s' + ^°''^) 

<^(^i + A;^2 + A;^°) + C'6V"*. (5.5) 
Using Gronwall's inequality again, we obtain from this that for < t < T 

E\t) + E\t) + 1 f e^'^^'-'^V\s)ds < Cl^^^K (5.6) 

Here we have used the fact that k is bounded and non-zero to absorb it into the constant C. 
We then employ (I4.55P with /3 = — 1 to see that 

=^ <S:!f (t) < £!f (0)e*^/2 + C f e^/^^'-'\v\s) + V\s))ds 

Jo 

^'-[(t) < Ct^e^^' £o"(i) < C£^-i{t) < Ci^e^^'. (5.7) 
Bootstrapping further, (I4.56P gives rise to 

E''{t) < Ct^e^^'. (5.8) 
Similarly, from f lCTj) . fl4T0|) . f l4TT]) . fl4T3|l and IKWf . we also obtain for < t < T 

£tit) + ^I'^it) + s'^'^'it) + £l::{t) + £l:;;{t) < c^^^\ (5.9) 

Next, from (I4.75P we get 

< (r/ + CQr)E^ + C^^^' £\t) < Ci^e^^K (5.10) 
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In turn, from f l4.83p and f l4.85p . we find that 

SX{t)+£t{t)<Ci'e''\ (5.11) 
Similarly, the energy inequalities ( I4.86p . fl4.88p . f l4.90p yield 

S\t) + £l{t) + £!f (t) < CL''e^^\ (5.12) 

and (Oip and fOS]) yield 

£\t)+£l'''{t)<Ci^e^'^\ (5.13) 
Successively, (S35]), ( lOHll . and (HSTl) imply 

£\t) + £j(t) + ^I'^t) < Ci'e'^*. (5.14) 

To get the bound of £^, we first note that £^ satisfies the following inequality from f l4.98p 
and flCT]) 

^[^' + F,] <{v + C9,){£' + F,) + C^\£\' + Ci'e'^' (5.15) 



for some constants Ci > and C > 0. We now define T* by 
T* := sup < t 



£:(s) < min 1^1,^1 forsG[0,t]|. (5.16) 

Let < t < min{T, T*}. Then by the Gronwall inequality, f l5.15p implies that 

£\t) + F,{t) < Ct^e^^' £\t) + £l^^{t) < Ci^e^^^ for < t < min{T, T*}. (5.17) 

Thus, combining all of the above analysis, we finally obtain 

£{t) < C^i^e^^^ for < t < min{T, T*} (5.18) 

for a constant C2 > independent of i. 

We now choose such that C^iQ^^'^ < min{6'i, 4^}- We consider the following two cases. 

(i) T{L,6i,) < min{T, T*}. In this case, the conclusion follows without any additional 
work. 

(ii) T{L,e^) > min{T,T*}. We claim that it must hold that T < T* < T{L,e^), in which 
case the conclusion directly follows. To prove the claim, we note that otherwise we would 
have T* <T < T{l, 6^). Letting t = T% from (EH]), we get 

£{T*) < C^i^e^^^* < Cai^e^^^' = 02(6^)^ by the definition of T{l, 6^), (5.19) 

but this is impossible due to our choice of 6^, since it would then contradict the definition of 
T*. Since we then find our desired estimate in both cases, this concludes the proof of the 
proposition. □ 
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5.2 Further nonlinear estimates 



As a preparation of the proof of our main theorem, we recall the Navier-Stokes-Poisson 
system (14 ■2p can be written in perturbed form as in fl3.38p and (13.391) in terms of a and 



w := r v: 



N2. 



with the boundary conditions 

(^^^ (0, t) = and a{M, t) = , B{w) = at x = M, 
where the boundary operator B{w) is defined by (I3.34p . Nq is given as 



Nr. 



5 + — ) Airad^w — As 



r / 



- 1 



w 



x=M 



and A^i and N2 read as follows: 
Ni = -47r(2po + cr)adxW 



No 



2w' 



a 



Po 



Ml - M2, 



where 



Ml 



X 

i4 1 ^6 



— / -2^y 

Jo Po 



ci r 1 fa 



ro Jo P* \Po 



C2 



a 



^0 V^'o Po 



by Taylor expansion 



cr 



and 



where p*/po ~ 1 is a bounded smooth function of — 

Po 



M2 = 16n' ( 5 + y ) {rtd.ipod^w) - r'd,{pd,w)} 

167r^ {5+^]{irt- r'')d.,ipod.,w) - r'd^iad^w)}. 



(5.20) 



(5.21) 



(5.22) 



(5.23) 



(5.24) 



(5.25) 



It is possible to estimate these nonlinearities in terms of the energy S given by (I4.17p . 
We present these estimates now. 



Lemma 5.2. For each t, 



(B{Ni,N2) <C\£\^ and \Nb\ + IdtN^l + Id^N^l < C\S\' 



(5.26) 



where € is defined in (13.460 . 
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Proof. The second inequality follows directly from Lemma 14.91 For the first inequality, we 
only provide the details for the highest-order nonlinear term M2 in N2. Lower order terms 
may be estimated similarly. Throughout the proof we will write 61 to denote the left side of 
estimate fl4.100p : Lemma [4.91 them implies that 61 < Cy/S. 
By rewriting M2 as 



Mo 



167r2(5+ f 



^1 -1 

r 



—r^da^ipd^w) + r^a 
P 



(5.27) 



a 



-r^d^ipd^w) - r'^d^ [ - ] pd^w 
P VP/ 



a 



it is easy to see that 



/ 



^-^dx < cel£ < C\£\\ 



(5.2^ 



Next, 



167r2(5 + f; 



r / 



31 ^")'+;o \P\^d^^pd^^)+r'd^(Pjyd^ 



Ttpr^r 



Po^ 
P 



w 



+ 



- 1 



Hence from the definition of the energies and from the estimates in the previous section, 

(5.30) 

On the other hand, dtM2 reads as 

dtM2 
167r2(5 + f 



4 - -^r'^vd:,{pd:,w) - r^d:^{dtad^w) - r^d-,{ad^dtw) + {r^ - r^)9^(po5x<9^^i') 



Thus 



-A-r^d^{pd^w) - —r'^d^ipd^w) - r^d^ ( — ) pdxW - r^-d^{pd^dtw) 
r p \ P J P 

r'^d^ ( - ] pdxdtw + (rg - r'^)—d^{pd.j,dtw) + (r^ - r^)4 ( — ) pd^dtw. (5.31) 
\PJ P \PJ 

\dtM2\^ 



I ^^dx < Cel{£l': + £l^^ + £l) + C{£' + £ln£l''' <C\£\' 



(5.32) 
□ 
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5.3 Data analysis 

In order to prove our nonlinear instability result, we want to use the linear growing mode so- 
lutions constructed in Theorem 12. II to construct small initial data for the nonlinear problem, 
written in the perturbation formulation f l5.20p . Small data in the perturbation formulation 
correspond to initial data for f ll.17p - fll.20p that are close to the stationary solutions p = po, 
V = 0, r = tq. Unfortunately, due to the regularity framework (given by S as in fl4.17p ) in 
which we have proved our nonlinear estimates, we cannot simply set the initial data for the 
nonlinear problem (15.201) to be a small constant times the linear growing modes. The reason 
for this is that the initial data for the nonlinear problem must satisfy certain nonlinear com- 
patibility conditions in order for us to guarantee local existence in the energy space defined 
by S. Until now we have taken the local well-posedness theory for the nonlinear problem for 
granted, but we must now say a few words about the compatibility conditions in order to 
construct our desired initial data. 

Recall that we can rewrite the nonlinear problem fll.17p - fll.20p in the form fl5.20p - fl5.2ip 
with nonlinearities given by fl5.22p - fl5.23p . Let us concisely rewrite fl5.20p as 



dtX + CX = m{X) for X = 




(5.33) 



where 9T(X) is the nonlinearity given in terms of A^^i and A^2 by the right side of (15.200 . We 
will also rewrite the boundary conditions fl5.2ip as 

f{w/rl)U=o\ /^(ro^-r-^)U=o\ 
C:(X) := (tU=m = := OIb(X). (5.34) 

\B{w)U=mJ \ Ns J 

Here r is determined as a nonlinear function of a as usual. 

Rewriting the nonlinear problem as fl5.33p - fl5.34p now allows us to easily describe the 
compatibility conditions for the initial data. Given X(0) as initial data for X at t = 0, we 
can use fl5.33p to iteratively solve for c^X(O) for j > 1: 

dtX{0) = -CX{0) +m{X{0)) 

dfX{0) = -CdtX{0) + Dm{X{0)) ■ dtX{0) (5.35) 

= -c{-cx{o) + m{x{o))) + Dm{x{o)) ■ (-£X(o) + oi(x(o))), 

and so on for higher derivatives, where here D is the derivative of the nonlinearity. We may 
similarly compute c^9Te(X)(0): 

dmBiXm = D^siXm ■ dtX{0) = DmBiXiO)) ■ [-£X(0) + ai(X(0))]U=M, (5.36) 

continuing as above for higher derivatives. This procedure may be carried out indefinitely 
as long as X(0) is sufficiently smooth. However, we may also differentiate the boundary 
condition (15.340 with respect to time and then set t = to see that the data must satisfy 
the boundary conditions 

m^m = dimsixm for j > o. (5.37) 
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Since the terms c^X(O) and 9f0^g(X)(0) constructed in fl5.35l) - fl5.36l) are determined entirely 
by X(0), we then find that the data X(0) must satisfy the nonhnear compatibihty conditions 
given by substituting fl5.35p - fl5.36p into f l5.37p . 

For completely smooth solutions to the nonlinear problem, the compatibility conditions 
would have to hold for all j > 0. In our case, we only require solutions to remain in the 
energy space defined by £, and as such we must only solve for (9^X(0) for j = 1, 2, 3, given 
X(0). This then requires the compatibility condition from fl5.37p only for < j < 3. Of 
course, in order to guarantee that £^(0) is finite, we must have that (9^X(0), < j < 3, satisfy 
the integrability conditions in the definition of S{0). This in turn gives us a natural Hilbert 
function space IH with the following three properties. First, if X(0) G H, then we have the 
trace estimates needed to make sense of the boundary conditions in fl5.37p for < j < 3. 
Second, if ||X(0)||jj is sufficiently small, then £{0) < C||X(0)||g for some C > 0. Here 
the smallness assumption is needed to deal with the nonlinearities in fl5.35p - fl5.36p and the 
r terms in S. Third, the linear growing modes produced in Theorem 12.11 are in H. It is 
straightforward to extract the proper definition of H from S and to work out the details of 
the estimate of S{0); as such, for the sake of brevity we omit these. With HI defined in this 
way, it is then easy to use estimate fl2.13p of Theorem 12.11 in conjunction with fl2.8p - fl2.10p 
to see that the growing modes are in H. 

Now that we have stated the nonlinear compatibility conditions, we see why we cannot 
simply set X(0) = lXq with 

X„ = (^-J (5.38) 

for cr^ and = w^^/vq the growing mode solution constructed in Theorem 12.11 and l > a, 
small parameter. Indeed, these solve 



AXo + CXo = and C(Xo) = ^ ^{C^^o) = for all j > 0, (5.39) 

which in particular means that X(0) = lXq does not satisfy the nonlinear compatibility 
condition f l5.37p for j > 1. 

To get around this obstacle, we will use the implicit function theorem to produce a curve 
of initial data satisfying the compatibility conditions, close to the linear growing modes. To 
this end, let us define the map F : EI — t- M}"^ via 



F(X) 



[ ^b(X) \ 



(5.40) 



where we understand that c^X and (9^DTb(X) for j = 1, 2, 3 are computed in terms of X as in 
fl5.35p - fl5.36p . Let Xq be the linear growing modes as above and let Xj G H, i = 1, . . . , 12, be 
arbitrary for now, with exact values to be chosen later. We then define / : IR-'^^-'^^ — t- IR-*^^ via 



12 



/(t, r) = F tXo + TiXi for t G M and r G 



[)12 



(5.41) 



i=l 
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It is straightforward to see, given the structure of the nonhnearities and 



^^2). Also, /(0,0) = and 
dt 



(0,0) 



/ ^(:^o) \ 

c:(AXo) 

^(A^Xo 



df 

and 7T^(0,0) 



that 



(5.42) 



From this it is then straightforward to choose the j£j for i = 1, ... ,12 so that the 12 x 
12 matrix 9//9r(0,0) is invertible. The imphcit function theorem then provides a small 
constant 6o > and a function ^ : (—to, '•o) 1^^^ so that f(t,C,(t)) = for all t G (— '•o) 
and so that ^ G and ^(0) = 0. We may then differentiate the equation f(t,C,(t)) = with 
respect to t, set t = 0, and use the first equation in fl5.42p to see that 



= ^(0.0) + ^(0.0)^ 







(5.43) 



since the matrix df/dT{0, 0) is invertible. Then ^ G with ^(0) = ^(0) = so that ^(t)/t^ 
is well-defined and continuous on (— ^o, '•o)- Using this, we may then deduce the existence of 
a small parameter Lq > and a curve ^ : {—lq, lq) ^ M given by 



12 



(5.44) 



i=l 



SO that for all l G (— ^o, '■o) 



F(2)(6)) = 0, i.e. 2)(t) satisfies the nonlinear compatibility conditions 

^/£m))<Cmi)\\^<CL, and (5.45) 

C($^(6)i,?)(6)2) <C, 

where here the norm ||-||q < is given by f l3.45p . the term (E is defined by fl3.46p . and in 
the second line we have written £^(2)(i)) for ^^(0) computed from the initial data X(0) = 2}(i). 
We now recast the above discussion as a lemma. 



Lemma 5.3. Let cr^,,Vi, be the growing mode solution constructed in Theorem \2.1l write and 



Wi, = r^Vi,, and assume the normalization 



(5.46) 



for ll'llg the norm defined by fl3.45p . Then there exists a number lq > and a family of initial 
data 



V^(0) 

for L G [0, Lq) so that the following hold 



x(0 



o-o(0 



(5.47) 



1. satisfies the nonlinear compatibility conditions required for a solution to the non- 

linear problem fl5.33p to exist in the energy space defined by £. 
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2. If £{0) denotes the value of S determined att = from the data then £{0) <Cl^ 
for a constant C > 0. 



3. For all l G [0, Lq), it holds that 



WqU) 



< £((To(0,w^o(0) < C 



(5.48) 



for a constant C > independent of l, where S is given by (I3.46P . 



4- Let ip'' denote the function given by (13.401) with Nq = A^g(j£(i)) determined by the data 
X(t) at t = 0. Then w''{0) = w''{0) — if)'' satisfies the homogeneous boundary condition 
B{w'{<S)) = and 

'0 



for a constant C > independent of l. 



(5.49) 



Proof. Everything except for the last item is proved above. The last item follows from 
Lemma 13.41 and the fact that Nq is at least a quadratic nonlinearity. □ 



5.4 Instability 

We are now ready to prove our main result. 

Theorem 5.4. There exist 6q > 0, C > 0, and < lq < 6q such that for any < l < lq, 
there exists a family of solutions <y''{t) and v''{t) to the Navier-Stokes-Poisson system (14. 2 p 
so that 

V£{0) < Cl but sup v^(t) > Oq. (5.50) 

0<i<T' 

Here T' is given by T^ = j\n ^ . 

Proof. We divide the proof into steps. At several points in the proof we will restrict the size 
of 6. Whenever we do so, we assume that l is also restricted so that < t < lq < 6. We will 
choose the value of 6q in the final step of the proof. 
Step 1 - Data and the solutions 

Let us assume that Lq is as small as the Lq appearing in Lemma 15.31 and then let X{l) 
for L < Lq he the family of initial data for the nonlinear problem (I5.33p - (l5.34p given in the 
lemma. For < l < lq, we now let (^,) be solutions to the Navier-Stokes-Poisson system 
( 15.33p -( !5.34p with a family of initial data 



t=o 



"^(°)^-X(0 = .f"^V^'H'!)- (5-51) 



w''{0) J V^*/ \Wo{t^ 



The solution satisfies \/£{0) < Cl. 
Note that since 



dy 



Po{y) + La^{y) + L'^ao{L){y) 



1/3 

(5.52) 
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a Taylor expansion and item 2 of Lemma 15.31 allow us to estimate 



ro(x) 



r'-(x, 0) 



L°° 



V 

+ 2 



r'-{x, 0) 



(5.53) 



L2 



for a constant > 0, independent of t. From this, the normalization (15. 46 p . and the 
estimate (I5.48p . we may assume that t < lq with lq small enough so that 



< v/^o,<x' (o) + £0,«'(o) + £i'^'(0) + ^/£0''■'{0) < 2i. 



(5.54) 



Throughout the rest of the proof we will let £{t) denote the total energy, defined by 
(I4.17p . associated to the solutions a' and at time t. 
Step 2 - Control of the energy 
Let us define the constant 



O I ^'^ 11. ||l/2 



(5.55) 



It will be useful in determining the time-scale in which instability begins. Indeed, we define 
Tby 

T := sup |s ^8^'"'{t)+£^^^'{t)+£^^^'{t) + ^fS^^ < (4 + 5o)6e^* for < t < s| . 

(5.56) 

The estimate (I5.54p guarantees that T > 0. Then by Proposition 15.11 and (I5.54p . there exist 
and 6'^ > such that for < t < min{T, T(t, 6'^)} (with T(i, 9^) given in the Proposition), 



te 



xt 



(5.57) 



Let us assume that 6 < 0^, which means that T'' < T{l, 6-^), and hence that the estimate 
(I5.57P also holds for < t < min{T, T'}. Let us further assume that 6 is small enough so 
that y/£{t) < C^e is small enough so that the right side of the estimate in Lemma 14.91 is 
smaller than 1/2. In particular, this implies that 



Po 



L°° 



1 

< - 

- 2 



(5.5^ 



for all < t < min{T, T^}. By further restricting 6 to decrease the bound of the terms in 
(I5.58p . and using the identities in fl4.2UI) and fl4.2ip . we can also bound 



(f°''^' (t)+£°''''(t)+^^'"'(t)) < 



< 2 (^°'"' (t) + {t) + E^'""' (t)) (5.59) 



for < t < = min{T, T'}. 

Step 3 - Linear estimates for a'' and w'' 

Notice that because of the estimate fl5.58p . the boundary condition w''/{r''y\x=o is equiv- 
alent to w'/rgla-^o- We can then modify the problem fl5.33p - fl5.34p to have the form fl3.4ip - 
fl3.42p , the latter of which has the homogeneous boundary conditions f l3.42p . This leads us 
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to consider ^'(0) = ^'(0) —ip'' as in Lemma [5. 3 [ which satisfies B{w''{0)) = at x = M. We 
then have that 



AC 



AC 







Then the solutions (^t) to (15.331) can be written as given in (13.440 : 



(5.60) 



M^)\_^tcfo\ 1/ 



Jo 



{t-s)C 



2(3) J 6 Jo \dtN^{s)rl/3j 



Here the nonhnear terms N^, N^, and N2 are defined in terms of w'' and a'' via (I5.22p and 

Theorem 13.51 together with the nonhnear estimates of Lemma 15.21 imply that if t < 
min{T, T(t, 6'^)}, then 



w''{t) J \w-, 



^2^tC f O-Ql'-j 



<C\£{t)f + C / e^^'''^S{s)ds 
Jo 

Jo 



(5.62) 



for a constant ^2 > independent of i. On the other hand, because of the estimates (13.471) 
and (I5.48p - (l5.49p . we may estimate 



^2^tc [ cro{L) 
Wo{l] 



AC 



< L^Ce^'^/e{ao{i),WoiL)) + Ce^V^(0,^0 < ^^^36^* (5.63) 



for a constant A3 > independent of l. Then we may then deduce from (15.620 . (15.630 . and 
the normalization (15.46^ that 



< Le^^ + L^A-^e^^ + A2{ie 



Xt\2 



(5.64) 



Step 4 - Control of the r energy 

We now turn to control of the term First note that 



d r u 
dt 2 



ro 



dx < 



1 



To 



2 \ 1/2 
dx 



w 



2 \ 1/2 
dx 



(5.65) 



which together with (I5.58P implies that 

d 



dt^ ^ ^ - I6V2 



w 



2 \ 1/2 
dx 



(5.66) 
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We may then argue as in fl4.2ip to see that 



V 

Po 



w 



dx < 



m2 , 4e_6 



5po Id.^w'l + yPo^o 



4 

ax < - 
- 9 



. (5.67) 



Combining fl5.66p and (15.671) with (15.641) . we find that 



d 



27 



dt^ ^ ' - 8V2 



IPoIIl' 



1/2 



< 



27 
8^2 



Po\\]!^ (^e^* + AsL^e^' + A2{Le^'f) (5.68) 



for < t < min{T, T*^}. Integrating this from to t < min{T, T*"} and employing (15.541) 
then yields the estimate 



I6V2X 

= Boite^') + {LA^){ie^') + A.ite^')^ (5.69) 

for < t < min{T, T^}, where Bq is the constant defined above in (I5.55p . and ^4,^5 are 
constants independent of t. 

Step 5 - The bound T' <T 

We now claim that if 9 is taken to be small enough, then T' = ^ In ^ < T. Suppose by 
way of contradiction that T' > T. Then the first bound in (15. 59 p . (I5.64p . and (I5.69P imply 
that 

^SO'-'it)+SO'-'it)+S^'-'it) + < (2 + 5o + LA^)iie^') + {2A2 + A^){ie'''Y 

< [2 + 5o + lA^ + (2^2 + A^)e]{Le^') < [3 + Bo](^e^*) (5.70) 

for t < T\ if we assume that 6 is small enough so that 6{2A2 + A^) < 1/2 and LQA4 < 1/2. 
For this choice of 6, we then find from the definition of T that T > T\ a contradiction. 
Hence T'' <T for ^ sufficiently small. 

Step 6 - Conclusion: instability 

We now define the part of the norm by 



00 



a 

Po 



1 

dx + - 



w 



dx. 



(5.71) 



Note that ||-||oo < IMIo that by the normalization (I5.46p . we have that the data satisfy 



:=CooG(0,l). 



(5.72) 



00 



Also, we may argue as in the derivation of (I5.59P to see that 

1 



V^O'-'(t) +^0,.' (^) > 



V2 



(5.73) 



00 



for < t < = min{r, T'}. 
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Let us now further assume that 9 is small enough so that < Coo/4 and 60^3 ^ ^00/4. 
We can then combine fl5T2D . ([5S2D, and flHIB^ to deduce that 



1 

"72 



00 
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AT' 



1 

> —=Le 
~ V2 



XT' 



o"o(0 



00 



00 



Wo{l) 



2 T'£ 

i e 



1 


^T'C 






00 V2 






00 



f^o(0 



V2 



1 

7! 



— Le 



XT' ( 



_ ,2 T'C ( "^OUJ 



+ e 



V2 



= i= (^Coo - e,A, - A,9') > ^9. (5.74) 

Setting 9q = {9Coo) /{2\/2), we find that (15.501) holds. This completes the proof of the 
theorem. □ 
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